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Auger Electrons Resulting from K-Capture 


FRANKLIN MILLER, JR. 
Rutgers University, New Brunswick, New Jersey 


(Received May 10, 1945) 


The most intense group of Auger electrons due to x-rays associated with K-capture in Zn® 
would have energies lying in a narrow band within about 300 ev of a maximum value of 7056 ev 
(#1p284). These would arise from K — L* transitions in the resulting Cu atom. Electrons from the 
K-—LM transitions would be about } as numerous, lying within about 300 ev of a maximum 
value of 8016 ev (Hp302). The soft electrons which are emitted copiously from Zn*® were found 
to have an absorption limit in collodion of 0.15+.03 mg/cm*, approximately 7500 ev. When 
samples were placed 0.6 to 1.3 cm from the screen of a screen wall-counter filled with Hz at 2 
cm Hg pressure, the upper limit (by inspection) was found to be at Hp302 +30 (8000+ 1700 ev). 
Scattering in the counter gas and absorption in the sample are believed negligible. These data 
are consistent with the hypothesis that the soft electrons from Zn* are Auger electrons resulting 
from K-capture. Similar results, in less detail, were obtained in deflection experiments with Fe. 


HEN an unstable isotope decays by K- 
capture, the resulting vacancy in the 
K-shell may be filled by an L, M, --- electron, 
and the characteristic x-rays of the product atom 
are emitted.' Indeed, one strong proof for K-cap- 
ture has been identification of the atomic number 


_of the atom emitting the x-rays, since if emitted 


as a result of K-capture the energy levels of the 
product atom are involved. It is not necessary, 
however, that all K-capture processes result in 
emission of an x-ray quantum. By the Auger 
effect, a radiationless transfer may take place, the 
single vacancy in the K-shell being replaced by 
two vacancies in the L-shell (K —L? transition), 
or by a vacancy in the L shell and one in the M 
shell (K—LM transition). The energies of the 
Auger electrons emitted in such transitions de- 
pend entirely upon the atomic number of the 
product atom resulting from K-capture, 


'L. W. Alvarez, Phys. Rev. 54, 486 (1938). 
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EXPERIMENTS WITH Zn® 


The 250-day isotope Zn®* decays by K-capture,! 
as well as by emission of positrons of maximum 
energy”? 0.4 to 0.7 Mev. A y-ray of 1.0 Mev is 
present, but the absence of many fast electrons 
shows it to be only slightly converted,‘ a result to 
be expected from the low value of Z*. Compton 
recoils from the y-ray are also observed.’ Ap- 
parently K-capture is favored as an alternative 
process to positron emission, since there is an 
anomalously high y/§-ratio, and annihilation 
radiation is present to an inappreciable extent.‘ 
Most of the Cu nuclei must be formed in an 
excited state, giving rise to the y-ray. The Cu Ka 


2 L. DuBridge et al., quoted in J. J. Livingood and G. T. 
Seaborg, Rev. Mod. Phys. 12, 30 (1940). 

3S. W. Barnes and G. Valley, Phys. Rev. 53, 946(A), 
(1938). 

‘J. J. Livingood and G. T. Seaborg, Phys. Rev. 55, 457 
(1939). 

5 R. Sagane, Phys. Rev. 55, 31 (1939). 
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x-radiation has been found! by means of differ- 
ential absorption measurements. 

Calculated Auger Energies 


The energies of the Auger electrons from the 
resulting Cu atom were calculated to the required 











COUNTS 
MIN 
800-4 2N65 
ABSORPTION iN 
600 + COLLODION 
400- 
ee i a 
° 
T T T T 
: MG 
05 10 AS .20 cut 
05 10 LS 2.0 FRINGES 


Fic. 1. Absorption in collodion of soft electrons from Zn®, 


degree of accuracy by the method of Ference,‘® 
using Slater's’? screening constant theory. Term 
values were taken from Siegbahn.* The most ener- 
getic K —L? Auger electrons have an energy given 
approximately by Ex—2E i =8954—2(933) 
= 7088 ev. When account is taken of the extra 
work required to remove the second Lin electron 
due to changes in screening, a correction factor 
AE =2E,—E,—E;=32 ev must be subtracted, 
leaving 7056 ev (/7p284 gauss-cm). If we assume 
extreme jj coupling, we find the lower limit for 
K—-L? will be Ex —2E1, —AE =6730 ev (Hp277), 
but judging from Ference’s experimental work 
with germanium, these are less numerous. Devia- 
tions from jj coupling will determine the detailed 
structure of the Auger levels, but will not ap- 
preciably affect the extreme values here listed. 
Thus we conclude that the K —L’ electrons are a 
nearly homogeneous group, lying within about 
300 ev of'a maximum value of 7056 ev. 
Similarly, the K—LM group, which is ex- 
pected® to be about 4} as numerous as the K —L’ 
group, lies within about 300 ev of a maximum 
value of Ex- (Et + E my. v) = 8016 ev (#1p302). 
The screening correction in this case is negligible. 
Direct measurements by Robinson and Cassie® 


6 M. Ference, Phys. Rev. 51, 720 (1937). 

7]. C. Slater, Phys. Rev. 36, 57 (1930). 

8M. Siegbahn, Spectroskopie der Réntgenstrahlen (Ver- 
lagsbuchhandlung Julius Springer, Berlin, 1931). 

®H. Robinson and A. M. Cassie, Proc. Roy. Soc. A113, 
282 (1928). 
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on the Auger electrons resulting from photo. 
ionization of Cu by Mo Kaé@ gave lines between 
11p286 and [/p292, consistent with the theoretical] 
values derived above. 


Apparatus 


A screen wall-counter'® was constructed having 
a counter wall 1.9 cm in diameter and 10 cm 
long. The wall was of copper or galvanized iron 
window screen, the insulation was of Lucite, and 
the wire was No. 30 gauge bare copper. The 
sample could be placed 1.3 cm from the screen, 
or closer if desired. Decay constants were checked 
by means of a thin glass-walled argon-oxygen 
counter. A Neher-Harper" quenching circuit and 
a scale-of-16 scaling circuit’ were used to record 
pulses. A cathode-ray oscilloscope was always 
used to check the behavior of the screen wall 
counter, which had a cosmic-ray background of 
about 40 per minute. 

Many very soft electrons were found to be 
emitted by invisibly thin zinc samples electro- 
plated on copper or aluminum foils. A sample 
giving about one count per minute (above a 
background of 8) on the glass counter gave about 
500 per minute in the screen wall-counter. The 
Zn® was over 2 years old at the time of the ex- 
periments (1943), and the decay curve, obtained 
from y's registered on the glass counter, was 
straight over two half-lives, showing a slope of 
258+10 days. 


Absorption Experiments 


Thin films of collodion placed 1 mm from the 
source were used in absorption experiments in the 
screen wall-counter. The chief error lay in de- 
termination of the thicknesses of the films, which 
were estimated in two ways: (1) from the 


TABLE I. Absorption in collodion of soft electrons from Zn®, 








Method Cut-off 
Color thickness 0.13 mg/cm? 
Fringes 0.15 
Weighing 0.2 








10 W. F. Libby and D. D. Lee, Phys. Rev. 55, 245 (1939). 

"1H, V. Neher and W. W. Harper, Phys. Rev. 49, 940 
(1936). 

2H. Lifschutz and J. L. Lawson, Rev. Sci. Inst. 9, 83 
(1938). 
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interference colors and (2) from the shift of 
fringes of Na light in a Michelson interferometer. 
Reasonable agreement was obtained. An_in- 
dex of refraction of 1.5 and a specific gravity 
of 1.66 were taken for collodion. A very rough 
check was made with a microbalance. The shape 
of the absorption curve is shown in Fig. 1; the 
points at 1.0, 1.5, and 2.0 fringes were determined 
by the interferometer and that at 1.0 fringe was 
also determined from the film’s color so as to 
provide overlapping of the two methods. Cut-off 
was at about 0.15+.03 mg/cm? (see Table 1). 
Extrapolating Schonland’s'* data for homogene- 
ous cathode rays, an Hp value of 285 was found 
(7500 ev). 


Magnetic Deflection Experiments 


The screen wall-counter was placed in a longi- 
tudinal magnetic field formed by an electromagnet 
from which pole pieces of diameter 10 cm had 
been removed. The coils were each 9 cm long and 
34 cm in outside diameter, and their centers were 
20cm apart. A rough approximation to Helmholtz 
coils was thus obtained. The counter was placed 
at the center of the system, and the magnetic 
field was explored by a search coil and found to be 
uniform to within 0.2 percent perpendicular to 
the axis, and to within 3 percent along the length 
of the counter. At cut-off a much shorter portion 
of the axis was used, indicating a more uniform 
field at cut-off than the outside limits of 3 
percent. The constant of the magnet was de- 
termined as 99+3 gauss/amp. with the aid of a 
search coil and a ballistic galvanometer. As the 
magnet current was increased, eventually only 
those electrons were counted which left the source 
in a direction perpendicular both to the field and 
to the line joining the source to the counter. Thus 
the radius of curvature was half the closest dis- 
tance of source to screen. A typical curve is 
shown in Fig. 2. It will be noted, as found by 
Libby and Lee,'® that a sharper cut-off was ob- 
tained with H. gas at 2 cm of Hg than with argon 
at 1 cm, because of scattering in the counter gas. 
In each case a trace of alcohol vapor was added 
to improve the counting characteristics. It is 
believed that errors due to scattering in the gas 


%O. D. Chwolson, Traité de Physique (Libraire Sci- 


entifique A. Hermann, Paris, 1906), Vol. 2, p. 596. 
* B, F. J. Schonland, Proc, Roy. Soc, A108, 187 (1925). 


and self-absorption in the source were negligibly 
small. The weighted mean of four runs, made at 
values of p from 0.27 to 0.64 cm, gave an end 
point (by inspection) of 302+30 gauss-cm, or 
Emax =8000+1700 ev. The rather large error 
arises from the difficulty of estimating end points 
(+5 percent), lack of homogeneity in the mag- 
netic field (+3 percent), and uncertainty in p 
(+5 percent). The sign of the particles was 
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Fic. 2. Magnetic deflection of soft electrons from Zn®, 
using screen wall-counter. A sharper cut-off is obtained 
when Hz: is used as the counter gas. Calculated Auger 
electron energies are indicated, those arising from the 
K—LM transitions being less numerous than those arising 
from the K —L? transitions. 


determined as negative, using sloping shields in 
the counter. 

From the shape of Curve B in Fig. 2, it would 
seem that the screen wall-counter can be used 
successfully at least down to 8000 ev electrons, 
which is considerably lower than the lowest 
energies measured by Libby and Lee. 


Discussion of Results 


Within the experimental errors, it is seen that 
both the absorption and the deflection experi- 
ments are consistent with the interpretation of 
these electrons as Auger electrons resulting from 
K-capture (See Table I1). 

Very few of the many experiments involving 
K-capture have been performed under circum- 
stances allowing detection of the extremely low 
energy Auger electrons. Plesset'® looked for 


' E. H, Plesset, Phys. Rev, 62, 181 (1942), 





TABLE II. Maximum Auger energies for Zn®. 











Method Maximum energy (ev) 
Theory 7056 (K—L?) 
8016 (K-—LM) 
Absorption 7500+ 1500 
Deflection 8150+1700 








electrons from Zn® in the range 20 to 300 kev and 
found none. Valley and McCreary"® searched for 
Auger electrons from K-capture in the 6.4 hr. 
Cd?!°7.1° and concluded that electrons of energies 
<20 kev have little action on the photographic 
plate. Hemmendinger"’ looked for x-rays from 
Mn®. He found a soft radiation which he was 
unable to identify, apparently because of absorp- 
tion in air. Unless the source was placed inside his 
ionization chamber, this radiation would be con- 
siderably harder than the Auger electrons. Hurst 
and Pool'® have found ‘electrons of very low 
energy associated with the fluorescence yield”’ in 
the K-capture of Ag'®*. No energy measurements 
were made, however. 

Auger electrons would arise following emission 
of any strongly converted nuclear y-ray. The 
justification for describing the Auger electrons of 
this paper as “‘resulting from K-capture”’ lies in 
the very small*'® internal conversion of the 


6G. E. Valley and R. L. McCreary, Phys. Rev. 56, 863 
(1939). 

17 A, Hemmendinger, Phys. Rev. 58, 929 (1940). 

LL. K. Hurst and M. L. Pool, Phys. Rev. 65, 60(A) 
(1944). 
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strong’ 1.0 Mev y-ray of Zn®. Practically the 
only way in which an electron leaves the K shell 
of Zn® is through capture. 


EXPERIMENTS WITH Fe* 


A small amount of Fe was available. This 
sample had been precipitated as Fe3;PQO,, and 
small traces of Mn and Co might have been 
present. At the time of the measurements, the 
sample was at least 4 years old, and showed a 
slow decay period of several years. Many soft 
electrons from this material were found in the 
screen counter, presumably from K-capture in 
Fe*®® (4 year half-life). A sample of purified Fe 
only a few months old also gave numerous soft 
electrons when a thin source was electroplated on 
a copper foil. A main group of maximum 
Hp220+30 (4300+1500 ev) was found, and 
traces of another group up to 20,000 ev possibly 
were also present. These results are based upon 
only one run. The calculated Auger maximum 
energies, however, are 5179 ev (K —L*), and 5872 
ev (K—LM), which are not discordant with the 
measurements. 
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Radioactive Scandium. II 


Car_ T. Hippon* anp M. L. Poor 
Mendenhall Laboratory of Physics, The Ohio State University, Columbus, Ohio 


(Received March 12, 1945) 


Sc?; This radioactive isotope has been reported to have a half-life of 2.62 days and to emit 
1.1 Mev electrons. These observations are not confirmed. A new radioactive isotope has, 
however, been produced by bombarding calcium with alpha-particles and to some extent by 
bombarding calcium with deuteron and proton. It emits beta-rays of 0.46 Mev and has a 
half-life of 3.4 days. The assignment is made to Sc*’. 

Sc**; This isotope has been produced by the reactions Ca**(p, m), Ca**(d, 2m), Ti**(m, p), and 
V"(n, a). From the V®(m, a) reaction gamma-rays of 1.33 Mev and beta-rays of 0.57 Mev 
were observed and from the Ca**(p, m) reaction gamma-rays of 1.35 Mev and beta-rays of 
0.57 Mev were observed. Relative saturation intensities have been calculated for Sc“, Sc, 
Sc*?, and Sc** when produced by proton, deuteron, alpha-particle, and fast neutron bombard- 
ments. Sc** was produced free of Sc“ by the V""(m, a) and the Ca**(p, m) reactions. Electrons of 
0.57 Mev and gamma-rays of 1.33 Mev were observed. Sc** emits approximately 14 gamma- 
rays per electron. All reactions associated with the scandium region are schematically sum- 


15, 1945 





marized in a nuclear transmutation chart. 


Sc‘? 


RADIOACTIVE isotope of scandium has 

been reported which decays with a half-life 
of 63+2 hours.' Beta-rays of 1.1 Mev were ob- 
served. Gamma-rays were detected, but it was 
not certain that the 63-hour period emitted 
them. This period was tentatively assigned to 
Sc*. 

In the present investigation, although several 
attempts were made to produce this isotope, it 
was not observed. 

A new radioactive isotope, however, has been 
produced by bombardments of calcium with 
alpha-particles. Various target holders made of 
copper, platinum, or gold were used in order to 
eliminate possible target-holder impurities in 
the scandium fraction. 

Figure 1 shows the results obtained for a short 
bombardment (4 hours) and a long bombard- 
ment (2 days); the half-life period is 3.40 days. 

Electrons of maximum energy 0.46 Mev were 
measured by absorption in aluminum. One of the 
absorption curves is shown in Fig. 2. Gamma- 
rays, if present, were too weak to measure. 

The assignment of this period to Sc*’ has 
been made for the following reasons: This radio- 
active isotope is most easily produced by bom- 
barding calcium with alpha-particles. Since Ca“ 


* Now at Cornell University, Ithaca, New York. 
1H. Walke, Phys. Rev. 57, 163 (1940). 
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is the second most abundant stable isotope of 
calcium, it is probably the one from which the 
new isotope is produced. Neither alpha-particle 
bombardments of potassium nor fast neutron 
bombardments of scandium produced the 3.4-day 
period. Therefore, an assignment to Sc® or Sc“ 
is improbable. The 3.4-day period could not be 
expected at Sc® or any scandium isotope of 
lighter mass, since no positrons nor gamma-rays 
were observed. Sc** is eliminated by the fact 
that Sc(d, p) and Sc“(n, y) reactions failed to 
produce the isotope. A V®'(m, a) reaction pro- 
duced only one period, the 1.83-day period 
which is obviously assigned at Sc**. 

The 3.4-day period has been also produced, 
although weakly, by bombarding calcium with 
protons and with deuterons. 

The following reactions have, therefore, been 
established: 

Ca“+a—Sc"+p 


Ca**+d—Sc?+n 
Ca**+ p—Sc"? +7. 


Sc*8 


A half-life period of 1.83 days was first assigned 
to Sc** by Walke.? Later by bombarding titanium 
and vanadium metals with fast neutrons he was 
able to make absorption measurements to deter- 


2H. Walke, Phys. Rev. 52, 777 (1937). 
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Fic. 1. Sc**: Decay curve showing a half-life of 3.4 days for beta-rays. 


mine the energies of the beta-rays (electrons) and 
of the gamma-rays. Two groups of beta-rays, one 
of 0.5 Mev (90 percent) and another of 1.4 Mev 
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Fic. 2. Sct?: Aluminum absorption measurements of the 
beta-ray spectrum obtained from the Ca(a, p) reaction. 


The end-point is,0.46 Mev. 
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Fic. 3. Sc**: Decay curve showing a half-life of 1.81 
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Fic. 4. Sc*: Aluminum and lead absorption measurements showing a beta-ray end-point of 0.57 
Mev, and a gamma-ray energy of 1.35 Mev. 


(10 percent), were reported. A gamma-ray of 0.9 
Mev energy was also reported. 

In the present work, long bombardments with 
fast neutrons were made. In order to increase 
the amount of vanadium in the target material, 
vanadium pentoxide was pressed into small 
pellets under great pressure. In this way sufficient 
activity was obtained to follow the decay for ten 
half-lives from which the half-life was found to 
be 1.83 days for both gamma-rays and electrons. 

By lead absorption measurements, a gamma- 
ray of 1.33 Mev was observed, which is the same 
energy as was found for the 3.92-hour period of 
Sc#. An aluminum absorption measurement 
shows that only one group of electrons of 
maximum energy of 0.57 Mev is present. 

Having determined the characteristics of Sc*® 
by the V°(m, a) reaction, which produced Sc** 
free of Sc“, it was then found that Sc** is also 
produced free of Sct by bombarding calcium 
with 5 Mev protons as shown in Fig. 3. For both 
the gamma-rays and electrons a half-life of 
1.81+0.02 days was calculated. Absorption 
measurements with aluminum indicate electrons 
of the same energy as for V*'(n, a) reactions. By 


absorption in lead the energy of the gamma-rays 
was found to be 1.35 Mev. Absorption measure- 
ments for both bombardments are shown in 
Fig. 4. 

By the method of semicircular focusing of 
Compton recoils in a magnetic spectrograph 
Mandeville has likewise found gamma-rays 
which he attributed to decay of Sc**. His measure- 
ments indicate a quantum energy of 1.35 Mev.* 

Electrons of 0.64-Mev energy attributed to 
Sc*® have been reported by Smith.‘ These 
measurements were made with a magnetic beta- 
ray spectrometer on scandium fractions from 
calcium bombarded with deuterons. 

Fast neutron bombardments of titanium 
produce weak Sc** activities. 

In all bombardments, the ratio of the intensity 
of the ionization produced by the beta-rays to 
that by the gamma-rays is approximately 15, 
which corresponds to an emission of about 14 
gamma-rays per electron. 

K-electron capture is therefore strongly sug- 
gested in lieu of many excited titanium nuclear 

*C. E. Mandeville, Phys. Rev. 62, 555 (1942); 64, 147 


(1943). 
* Gail P. Smith, Phys. Rev. 61, 578 (1942). 
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Fic. 5. Schematic summary of relative saturation intensities of various scandium 
periods obtained by various reactions. 


levels.’ This suggestion was given support by 
several absorption measurements which were 
made with the inverted electrometer and the 
thinnest pieces of aluminum and beryllium ob- 
tainable. Characteristic x-rays associated with 
K-electron capture for this region are of the 
order of 3A. The experimental values obtained 
ranged between 2.4A and approximately 3A. 


RELATIVE INTENSITIES 


All of the radioactive periods and their relative 
intensities are schematically shown in Fig. 5. 
Solid lines represent typical decay curves for 
beta-rays and broken lines, gamma-rays. 

Relative intensities, as calculated for infinite 
bombardment, are indicated on the arbitrary 
basis of 1000 for the most intense half-life period 
produced. The ratio of the intensity of the beta- 
rays to the intensity of the gamma-rays is 
indicated for each half-life period. 

Two main periods were produced by alpha- 
particle bombardments of calcium. The 3.4-day 
period of Sc‘? and the 3.92-hour period of Sc* 
have saturation intensities of 1.11 and 1000, 
respectively. The latter is the most intense 


5 Ernest Pollard, Phys. Rev. 54, 411 (1938). 


period of scandium and may have an initial 
intensity 30,000 times background. 

The isomers of Sc, when produced by the 
K*"(a, m) reaction, are found to have approxi- 
mately equal saturation intensities; namely 900 
for the 3.92-hour period and 1000 for the 2.44-day 
period. 

Sc*® was the main isotope produced by 5 Mev 
proton bombardments of calcium. Sc* was also 


TABLE I. Summary of results. 











Energy of Energy of 
Active Observed Produced beta-rays gamma-rays 
_tectope periods by in Mev in Mev 
Sc® Not 
observed 
Sc* 3.92 hours Ca-a-p 1.13(+) 1.65 
Ca-p-n 
Ca-d-n 
Ca-d-2n 
Sc# 3.92 hours K-a-n 1.33(+) 1.33 
Sc-n-2n 
Ca-d-n 
Ca-d-2n 
2.44 days 0.28 
Sc‘? 3.40 days Ca-a-p 0.46( 
Ca-d-n 
Ca-p-y 
Sc** 1.83 days V-n-a 0.57(—) 1.33 
Ti-n-p 
Ca-d-2n 


Ca-p-n 
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Fic. 6. Transmutation chart showing the scandium region and reactions leading to radioactive scandium isotopes. 


produced, but no evidence was found to indicate 
Sc; however, a weak period of 3.4 days due to 
Sc‘? was present. The saturation intensities of 
these isotopes, Sc**, Sc, Sc, and Sc*’, are, 
respectively, 1000, 570, 0, and 80. 
Deuteron-bombardments of calcium produced 
several radioactive isotopes of scandium, namely 
Sc, Sc“, Sc*7, and Sc*’. Sc#8 was the main isotope 
produced and was approximately seven times as 
intense as Sc“. The composite decay curve of 
these activities is so complex that an attempt to 
discover the individual periods by subtractions 
could only fail. The half-life discrepancies re- 
ported in the literature based upon deuteron 
bombardments may in this way be reconciled. 
The saturation intensity of the Sc* 3.92-hour 
period plus the Sc“ 3.92-hour period is 700. The 





intensities of the 1.83-day period of Sc** plus the 
2.44 day period of Sc“ is 1000. The 3.4-day 
period of Sc*’ has an intensity of 83. 

The results of this investigation and that 
reported earlier® are summarized in Table | and 
Fig. 6. 
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A two crystal vacuum spectrometer was used to investigate the x-ray absorption spectrum of 
solid potassium chloride over a fifty-volt range in the vicinity of the K edges of both potassium 
and chlorine. Absorbers were prepared by vacuum evaporation of the salt onto thin, taut 
Kodapak films. High precision of intensity measurement was made possible by the use of 
fortuitously located La; radiations from targets of antimony and palladium. No corrections to 
observed intensities were made for the finite resolving power of the instrument. Results are 
expressed by graphs of relative absorption coefficient as a function of relative photon energy, 
and are summarized in a table of relative energy separations of prominent absorption maxima 


and minima. 





INTRODUCTION 


HE study of fine structure in the x-ray 
absorption edges of solids has long been 
considered important for an understanding of 
their electronic states. Recently, successful corre- 
lations between theory and experiment have been 
made for some of the simpler metals.'? However, 
both theoretical and quantitative experimental 
work on non-conductors is meager and unsatis- 
factory. The present experimental work, using 
the most precise techniques known, was under- 
taken in order to provide information about the 
x-ray absorption fine structure of potassium 
chloride, an ionic crystal which is presumably one 
of the simplest from a theoretical standpoint.’ 
The experimental equipment consisted of a 
high voltage d.c. power supply, an x-ray tube, a 
set of potassium chloride absorbing films, a two 
crystal vacuum x-ray monochromator, an ioniza- 
tion chamber for detecting the x-rays, and a d.c. 
amplifier with galvanometer output for measuring 
ion currents. Several runs of absorber opacity 
versus monochromator setting were made for 
each absorption edge. The natural logarithms of 
the opacities for each run were then suitably 
corrected and the values for separate runs ad- 
justed and matched so as to yield results con- 
veniently represented as graphs of relative 
absorption coefficients as a function of relative 
photon energy. 


1W. W. Beeman and H. Friedman, Phys. Rev. 56, 392 


(1939). 
2 J. B. Platt, Phys. Rev. To be published. ; 
*Salient features of the fine structure of potassium 
chloride have been reported by J. M. Nuttall, Phys. Rev. 
31, 742 (1928); and G. P. Brewington, Phys. Rev. 46, 861 


(1934). 


EXPERIMENTAL DETAILS 
Apparatus 


The high voltage power supply was of the full- 
wave rectifier type with filtered and electronically 
stabilized output.‘ The x-ray tube voltmeter, 
consisting of a galvanometer in series with two 
1000-megohm resistors could be read with a 
precision of +0.1 percent and with an accuracy 
of +0.5 percent at 10,000 volts. The tube current 
was read with a 0-100 milliammeter with a pre- 
cision of +0.2 percent. The majority of measure- 
ments were made with a tube voltage of 10,000 
volts and current of 100 ma. 

Details of the x-ray tube construction have 
been given elsewhere.’ The target used for the 
study of the potassium edge was an alloy of 12 
percent antimony and 88 percent lead tinned to 
the copper target carriage,’ while the target used 
for the study of the chlorine edge was a 10-mil 
sheet of palladium soldered to the carriage. 
These target materials were selected because the 
La, lines of antimony and palladium provided 
both high x-ray intensity in the regions of the 
edges and standards for wave-length meas- 
urement. 

The absorbing films were prepared by vacuum 
evaporation of chemically pure potassium chlo- 
ride from a furnace, similar to one described by 
Schneider and O’Bryan,® onto tautly stretched 
Kodapak films 0.001 inch thick. The thickness of 


‘L. G. Parratt and J. W. Trischka, Rev. Sci. Inst. 13, 17 


(1942). 

5 L. G. Parratt, Phys. Rev. 54, 99 (1938). . 

6 E. G. Schneider and H. M. O’Bryan, Phys. Rev. 51, 293 
(1937). 
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TABLE I. Energy separations in electron volts of fine 
structure maxima and minima from initial absorption 


peak, 











Feature K —K edge Cl—K edge 
A 0.0 0.0 
a 1.7 y 
B 3.9 3.8 
B 5.5 5.2 
Cc 6.9 5.7 
Y ? 8.2 
D 9.0 10.1 
6 11.5 11.4 
E 16.3 14.3 
€ 22.4 ? 
F 27.5 19.2 
i? — 23.4 
G — 29.9 








absorbers were made in the region of the potas- 
sium edge. Of these runs three covered the entire 
region studied. In the region of the chlorine edge 
three runs, with three different absorbers, were 
made. 

The general method of treating absorption 
measurements where more than one element is 
present in the absorber has been described by 
Platt.2 The jump ratios of both potassium and 
chlorine were taken to be 10. Successive runs 
were matched at the peak of the first absorption 
line. The corrections for higher order radiation, 
made to the natural logarithm of the opacities, 
had a maximum value of 3 percent for the 
potassium measurements and 1 percent in the 
region of the chlorine edge. No correction for the 
resolving power was attempted since no reliable 
method has been developed for making such 
corrections to absorption spectra. 

For purposes of absolute wave-length determi- 
nation, separate runs were made to determine the 
position of the first absorption maximum relative 
to the La; radiation used for recording the fine 


structure. 
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RESULTS 


Figure 1 is a graph of relative values of ux vs. 
relative photon energy in electron volts in the 
regions of the K absorption edges of chlorine and 
potassium in solid potassium chloride. The wave- 
length scale used in expressing the results is an 
X.U. scale rather than a Siegbahn scale. The 
standard wave-lengths used were 3431.8 Siegbahn 
units for Sb La; and 4358.5 Siegbahn units for 
Pd La,.* The ratio of the grating and Siegbahn 
scales of wave-lengths was taken as 1.002034," 
The wave-lengths of thie peaks of the first ab- 
sorption lines in the chlorine and potassium 
spectra were determined to be 4387.6 X.U. and 
3432.5 X.U. respectively." In establishing the 
voltage scale the wave-length associated with one 
absolute volt was taken to be 12,395 X10-* cm 
absolute volt.'® 

For ease of comparison with other data and 
with theoretical results, which are likely to take 
the form of energy separations, Table | presents 
the energy separations from the initial maximum 
of the various maxima and minima in Fig. 1. 
Maxima and minima are referred to by letters of 
the Roman and Greek alphabets respectively. 

It is estimated from the self and mutual con- 
sistencies of the various runs that variations in px 
as small as 3 percent are significant. An addition 
of component errors gives a maximum uncer- 
tainty in absolute wave-length specification of 
+0.15 X.U. Relative wave-lengths are precise to 
+0.02 X.U., the error in matching different runs. 

The author wishes to express his indebtedness 
to Professor L. G. Parratt for many stimulating 
discussions and helpful criticisms, and to Pro- 
fessor H. A. Bethe for suggesting the investigation. 


*H. Hjalmar, Zeits. f. Physik 3, 262 (1920). 

1 R. T. Birge, Rev. Mod. Phys. 13, 233 (1941). 

" Or 4378.7 Siegbahn units and 3425.5 Siegbahn units 
respectively. 
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By the theory of Sommerfeld, relativity effects and 
retardation of potential being neglected, matrix elements 
and associated components of continuous spectrum x-radia- 
tion are computed for a variety of electron-nucleus colli- 
sions and for a distributed series of positions in the spec- 
trum. The calculations cover values of V/Z? from 0.06128 
to 3.356, where Z is the atomic number of the nucleus, and 
V is the bombardment potential in electrostatic units. 
Accuracy of calculation is 1 percent. Screening is neglected 
except at the long wave limit of the spectra where it is 
taken into account by Sauter’s method. Empirical algebraic 
formulas are found which closely represent the rigorously 
calculated results. Intensity and polarization predictions 


for any direction of emission and any excitation conditions 
within the range of applicability of the theory may be 
readily drawn from the computed results. Elwert’s proposed 
correction factor for rectifying the approximate spectral 
intensities of Sommerfeld and Maue is found effective 
within the limits of its restricting assumptions. Theoretical 
efficiencies of continuous x-ray production are calculated 
by combining theoretical intensities with known rates of 
electron energy loss in traversing matter. Thick target 
efficiency (a ratio, not percent) is given by 1.4 10~*Z kv. 
Thin target efficiency is found to be approximately twice 
the thick target efficiency for any given Z and kv. 





HE problem of the continuous x-ray spec- 

trum is to find answers to such questions 
as the following: When an electron of stated 
kinetic energy collides with a thin layer of matter 
of atomic number Z and of some given surface 
density, what is the probability that a photon 
whose frequency lies between any stipulated 
limits shall be emitted in an elementary solid 
angle of specified magnitude oriented in any de- 
fined manner with respect to the direction of 
motion of the bombarding electron? A further 
question concerns the state of polarization of 
such a photon. 

Answers to these questions are implicit in 
extant theory, but in only a very few cases have 
explicit numerical answers been extracted from 
the theories and placed on the record. It is not 
yet possible to say broadly that the theory either 
does or does not agree with experimental facts, 
though a few encouraging similarities have ap- 
peared. We shall not have a satisfying number of 
such comparisons without further work on both 
sides of the gap, i.e., until many more good 
measurements have been completed and until 
theoretical predictions are made available in 
numerical form. 

The present paper attempts to supply the 
second of these necessities. It contributes nothing 
to the basic theoretical picture according to 
which the emission of the continuous x-ray spec- 
trum is now understood; rather, it presents a 
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computational procedure for obtaining numerical 
values of x-ray energy and polarization from the 
theory, together with the results of such cal- 
culations. 

For a system consisting of a bare nucleus and 
one electron, Schrédinger’s equation predicts a 
continuous range of positive energy states, and a 
continuous spectrum must result from the joint 
operation of the possible transitions between such 
states. The spectral features may be deduced, as 
Oppenheimer! and Sugiura? have shown, by 
calculating matrix elements from characteristic 
functions corresponding to states involved in the 
transitions, but the method introduced by 
Sommerfeld* and developed by a number of sub- 
sequent theorists has greater analytical simplicity 
and will be followed in the present paper. 

The theory of Sommerfeld considers the wave 
systems of an electron approaching an atomic 
nucleus and departing with altered direction and 
reduced speed, having suffered an energy loss 
which appears as the energy of the emitted 
photon. The single process can be described with 
rigor, but the integration over all possible direc- 


‘J. R. Oppenheimer, Zeits. f. Physik 55, 725 (1929). 

? Y. Sugiura, Sci. Pap. Inst. Phys. Chem. Res. Tokyo 17, 
89 (1931); and earlier contributions referred to therein. 

*A. Sommerfeld, Ann. d. Physik [5] 11, 257 (1931); 
O. Scherzer, Ann. d. Physik [5] 13, 137 (1932); F. Sauter, 
Ann. d. Physik [5] 18, 486 (1933), [5] 20, 404 (1934): 
A. Sommerfeld and A. W. Maue, Ann. d. Physik [5] 23, 
589 (1935); G. Elwert, Ann. d. Physik [5] 34, 178 (1939); 
R. Weinstock, Phys. Rev. 61, 585 (1942), 65, 1 (1944). 
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tions of departure of the electron—a necessary 
operation if theory is to be compared with ex- 
periment—has been found troublesome by all 
who have tried it. 

Approximations, physical or mathematical, 
have been necessary features of all the treat- 
ments. The present results are conditioned by the 
inexact physical assumptions of the original 
Sommerfeld* theory, but we follow Weinstock? in 
avoiding mathematical approximations such as 
those which restricted the applicability of the 
results of Sauter and of Elwert.‘ 

The questionable physical assumptions and 
approximations are the following: (a) Except in 
connection with the long wave limit of the spec- 
trum, atomic electrons are ignored, and the 
nucleus is assumed to be the center of a pure 
Coulomb field. (b) Relativity effects, retardation 
of potential, and the effects of electron spin are 
neglected. (c) The de Broglie waves of the ap- 
proaching and receding electron are treated as 
plane waves, filling all space. 

Picture a row of electrons moving along or 
closely parallel to the x-axis of a cartesian coordi- 
nate system (see Fig. 8), approaching the origin 
from the negative side with speed v; conferred by 
a potential difference V. Upon, arriving at the 
origin, the electrons strike or pass through a 
plane target coinciding with the yz-plane which is 
composed of atoms of atomic number Z in a 
uniform distribution. We investigate the con- 
tinuous x-radiation resulting from collisions at 
the origin and departing in a direction lying in 
the xz-plane and making an angle @ with the 
positive x-direction. The electron leaves the 
collision with speed v2. As applied to such a 
picture, the theory deals in probabilities of 
photon emission, but with large numbers of 
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atoms or of electrons, such probabilities increase 
proportionately, and one may speak with sta- 
tistical accuracy of emitted intensities, powers, 
and energies. 

The bombarded and emitting atoms simulate a 
source possessing calculable component Hertzian 
dipole moments, each of which produces radia- 
tion components J,, J,, and J,, proportional to 
the square of the corresponding moment. The 
total intensity proceeding in a direction @ con- 
sists in general of contributions from all three 
radiation components weighted with reference to 
the angles between their dipole axes and the 
direction of observation. The observable in- 
tensity is given by® 


I,=I, sin? 6+],+], cos? 6. (1) 


The calculation and discussion of J, ,, . the 
objects of this paper; these components will be 
defined in such a way that J» (not rigorously a 
radiation intensity) comes out in the units ergs 
per steradian per unit frequency range per bom- 
barding electron per atom-per-square-centimeter 
of target area. Neglect of electron spin leaves the 
direction of bombardment an axis of complete 
necessarily. There remain 


are 


symmetry, so I,=T, 
two quantities, 7, and J,, to calculate, and each 
is a function of the four parameters 6, V, Z, and 
the frequency v, a group which will be found 
reducible to three. In terms of these components, 
the intensity and polarization of radiation emitted 
in the direction @ are given by 


Ig=T, sin? 6+1,(1+ cos? 6) | 








I.y2= 


where N?, ,,, are dimensionless quantities equiva- 
lent to M?, ,,./|A|? of Weinstock’s first reference. 
The computation of these quantities is almost the 
whole of the problem confronting the computer. 


limited by Born’s first approximation, 


‘Sauter is 
while Elwert’s method is restricted by the 


Ze*/hy, Kil 4 


[1—exp (—42*e?Z/veh) |L[exp (49°e?Z vh) —1){(I ’e)\— (Ve—hyv) nc? 


1—(1,/Iz) +. (2) 
“14, [,)(2 csc? @—1)| 
The method of Sommerfeld yields: 
4mre®Z?( Ve—hv) N*,y,: 
oe (3) 
CALCULATION OF N,? (=N.’) 
From Weinstock’s derivation one finds 
16°Z* -e3m 29°Z*e'm 
ee, ae 
h? VE,° h? (Ve—hp) 
.. w* . 
xf | ae | |G3|"dw. (4) 
1 a (1—w)? 


assumption (Ze*/hve)—(Ze*/hn)<K1 or Zeé/hvy,<Ki. The 


effects of these conditions will be noted on a subsequent 
page. 


G. ~ 8G, Elwert. reference 3; p. 194. 
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New symbols appearing in Eq. (4) are defined as 
follows: 


fo= —4nyn2/(mi—M2)*, wo= ko/(Eo—1). 


The quantity G; is the hypergeometric function 
F(i+m, —me2, 2, w) where m,; and me are imagi- 
nary quantities defined by m,2=(—2me*Z1) / (hv, 2). 
The physical significance of w is of no concern 
since it goes out in the definite integration. 

Weinstock expanded Gz; in series and inte- 
grated Eq. (4) term by term, obtaining his 
compact Eq. (9). Although this works well in 
some cases, it has not been found generally 
trustworthy because of the danger that the series 
of integrals may be chopped off too soon. It is 
quite impossible to depend upon appearances in 
estimating the progress of its convergence, and 
no positive criterion has been found. It is known 
that as many as twenty terms are insufficient in 
certain cases to reduce the residue to negligibility, 
so the evaluation is both suspect and laborious. 
It was decided, therefore, to evaluate the inte- 
grand in Eq. (4) for representative values of w 
and to carry out the integration graphically. 

The range of integration extends from zero to 
Wo, a positive quantity lying between zero and 
unity. For those cases (values of V and v) which 
give values of w» not too near to unity, say not 
exceeding 0.75, the necessary values of G3 may be 
obtained from the standard hypergeometric 
expansion 

abx a(a+1)b(b+1)x? 
Pia, 6, ¢, s)a14-—4——______—_ 
lc 1-2-c(c+1) 
a(a+1)(a+2)b(b+1)(6+2)x* 


1-2-3-c(c+1)(c+2) 





u=x 
“= ® gux" 
u=0 


where » is an integer designating a series term, 
which starts with 4» =0 for the first term. The 
coefficients are given by 


go=1, and in general 
(u+a)(u+d) 

a Ru: 
(u+1)(u+c) 





Rutt 


By use of this series, G; may usually be evaluated 
with an uncertainty not greater than one percent 
without going beyond 10 or 15 terms. For values 


of w too large for prompt convergence, use is 
made of the identity® 


G3;=C,F(1+n,, —Neo, N1— Ne, 1—w) 
+C2F(1—,, 2+, 2—ni+n2, 1—w), 





where 
—nNy+Ne , 
C=———_—4, 
—nno(1+72) 
A 
C,= —_——_——-- (1—w)!oites, 


n;(—1 4+-%,— 1s) 
A= (m,—n2) T'(m,)T(—nz), 


and A designates the complex conjugate of A. 
Numerical evaluation of A is most readily 
carried out by means of the identity 


1/I'(x) =x(x+1) 
XK (1+Bix+Box?+Byx'+---). (5) 


The coefficients of x as calculated by Bourguet’ 
are reproduced in Table I. 

The terms in the series of Eq. (5) are alternately 
imaginary and real, leading to a complex A. For 
the largest values of m2, the tabulated coefficients 
are insufficient, and it is advantageous to use 
Legendre’s duplication formula 


w'T'(2x) = 22*-'P(x)P(x+4), 


to bring the arguments of the gamma-functions 
within the range of the table. The imaginary 
character of m, and mz requires special attention 
in the evaluation of C2; here we may write 


(1 —w)!-"!+"2= (1 —w)(cos +47 sin 8), 
in which 
6= |n,—ne\ log (1/1—w). 


lf the integrand of Eq. (4) be plotted as a 
function of w, the curve will usually be found 
inconvenient for planimeter integration because 
of the high and narrow peak at the upper limit. 
This feature is avoided by plotting as a function 


* E. T. Copson, An Introduction to the Theory of Functions 
of a Complex Variable (Oxford University Press, New York, 
1935), p. 251. 

7L. Bourguet, Acta Math. 2, 261 (1883). Harold T. 
Davis, Tables of the Higher Mathematical Functions (The 
Principia Press, Inc., Bloomington, Indiana, 1933), Vol. 1, 
p. 185. 
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of u, where u=1/(1—w), in which case the 
errors of area measurement may well be less than 
those of computation. The transformed integral 
is 





rr (w—1)(fo-+u—1) 
|Gs| %du. 


u2 


As & increases, the maximum of the integrand 
becomes higher and narrower, but plotting in 
two or three different sections with suitable 
changes of scale readily keeps it in hand. 


CALCULATION OF N,? 


The calculation of NV,? might be carried out in 
a manner analogous to that described above, but 
Weinstock’s expression contains three definite 
integrals, and although he has shown how to 
replace them by infinite series, the problem of 
convergence persuades us to adopt a differing 
approach. Sommerfeld and Maue have shown 
how to express the sum N?=N?Z+N,?+N2 in 
closed form, and Weinstock has put their ex- 
pression into a form based upon four rapidly 
converging series, two of which we find to be 
identical. We therefore employ Weinstock’s equa- 
tion for N? in the form shown below and later 


TABLE I. Coefficients in the power-series expansion 


ain =x(x+1)(1+Bix+Box?+---). 











I'(x) 

n Ba 
1 — 0.4227 8433 5098 4671 
2 — 0.2330 9373 6421 7867 
3 0.1910 9110 1387 6915 
4 — 0.0245 5249 0005 4000 
5 — 0.0176 4524 4550 1443 
6 0.0080 2327 3022 2673 
7 — 0.0008 0432 9775 6044 
8 — 0.0003 6083 7816 2548 
9 0.0001 4559 6142 1399 
10 — 0.0000 1754 5859 7517 
1] — 0.0000 0258 8995 0224 
12 0.0000 0133 8501 5466 
13 — 0.0000 0020 5474 3152 
14 — 0.0000 0000 0159 5268 
15 0.0000 0000 6275 6218 
16 — 0.0000 0000 1273 6143 
17 0.0000 0000 0092 3397 
18 0.0000 0000 0012 0028 
19 — 0.0000 0000 0004 2202 
20 0.0000 0000 0000 5240 
21 — 0.0000 0000 0000 0140 
22 — 0.0000 0000 0000 0067 











solve for N,2=N?—2N,?. This avoidance of a 
direct attack upon JN,’ follows, in principle, the 
example of Elwert. 

The working equation for N? is 














—167f|A\? « 
N= — rqrég*—' 
£0 m1| =! 
~}-—(- enn] Y rite 
Nino r=] 





« 

—(m,—m2) inl] (6) 
r=0 

wherein #=“‘real part of,” and the quantities gq, 

and j, are themselves defined by the finite series: 


; * f . * 
¢, = > £4, £4r-+5 Jr= > £4, x25, r—pv- 
v=0 yv=0 


Quantities in these series are defined and calcu- 
lated by the recursion relations: 








£4,0>=25,0= 1, 
(v—1—m,)? 

£4 >= £4, v—1)5 
v(v-+ne—n}) 
(v-1 — mn)? 

£5.0> £5, v—1- 
v(v-+n,—n2) 


The asterisk (*) denotes the conjugate function. 
In these series vy and r are real integers assuming 
consecutive values within the stated ranges. The 
convergence is such that as few as three terms 
sometimes suffice to fix the value of a series. The 
g's are, of course, complex, but the conjugates in 
gr and j, cancel out the imaginaries and leave 
these quantities and their series real. 


NUMERICAL EVALUATION 


Inspection of the foregoing equations shows 
that the variables V, Z, and v are present only in 
combinations reducible to V/Z? and v/v. This 
reduction of the number of parameters con- 
trolling the spectral features from four to three 
(the two quotients above and @) greatly reduces 
the amount of computation necessary to cover 
the full range of experimental variables, since any 
one complete calculation furnishes a separate 
prediction concerning each of the 92 elements. 
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TABLE II. Components and total energy of radiation. 











V/Z? v/ vo Iz X10" Iy( =I.) X10 W x10” 
3.356 0.403 0.0809 0.0296 1.17 
2.171 0.623 0.128 0.0228 1.45 
1.776 1.00 0.146 0.00132 1.24 

0.761 0.154 0.016 1.56 
0.197 0.144 0.105 2.97 
0.0962 0.136 0.151 3.67 
0.00 0.0586 0.223 4.23 
1.351 1.00 0.207 0.00223 1.77 
0.7896 1.00 0.398 0.00761 3.47 
0.624 0.377 0.072 4.37 
0.443 0.361 0.121 5.05 
0.217 0.332 0.227 6.59 
0.00 0.132 0.496 9.42 
0.4934 1.00 0.685 0.0191 6.06 
0.708 0.637 0.0948 6.92 
0.346 0.578 0.254 9.10 
0.00 0.211 0.776 14.8 
0.3493 1.00 0.999 ~° 0.0378 9.01 
0.2961 1.00 1.20 0.0518 10.9 
0.577 1.05 0.250 13.0 
0.00 0.351 1.26 24.1 
0.1709 1.00 2.12 | 0.141 20.1 
0.1359 1.00 2.63 0.214 25.6 
0.761 2.43 0.402 27.1 
0.197 1.91 1.39 39.3 
0.0962 1.79 1.90 46.8 
0.00 0.765 2.705 51.7 
0.06128 1.00 5.43 0.697 57.2 
0.761 5.18 1.05 61.0 
0.197 4.26 2.92 84.5 
0.0962 3.92 3.95 99.0 
111.0 


0.00 1.70 


5.74 


| 
| 
| 


V/Z? v/vo 12X10" J,( =I.) X10" W x10" 
0.06041 1.00 5.51 0.718 58.2 
0.624 5.12 1.35 65.4 

0.443 4.86 1.79 70.6 

0.217 4.40 2.81 83.9 

0.00 1.72 5.81 112.0 

0.03776 1.00 8.48 1.30 92.9 
0.708 8.11 1.92 100.0 

0.346 7.34 3.32 117.0 

0.00 2.78 9.10 176.0 

0.03546 1.00 9.00 1.40 98.9 
0.761 8.73 1.92 105.0 

0.197 7.35 4.76 141.0 

0.0962 6.75 6.34 163.0 

0.00 2.95 9.60 186.0 

0.02723 1.00 11.5 1.92 128.0 
0.624 10.8 3.01 141.0 

0.443 10.4 3.84 152.0 

0.217 9.74 - 5.68 177.0 

0.00 3.81 12.2 236.0 

0.02265 1.00 13.62 2.37 154.0 
0.577 12.73 3.81 171.0 

0.00 4.60 14.5 282.0 

0.01702 1.00 17.7 3.27 204.0 
0.708 17.2 4.31 216.0 

0.346 16.0 6.84 248.0 

0.00 6.09 18.8 366.0 

0.01576 1.00 19.0 3.58 219.0 
0.624 18.2 5.14 238.0 

0.443 17.7 6.36 255.0 

0.217 16.6 9.25 294.0 

0.00 6.58 20.2 394.0 

0.009854 0.708 28.9 7.48 368.0 
0.346 27.2 11.2 415.0 

0.00 10.5 30.5 599.0 


| 
| 
| 
| 
| 





Many of these predictions, however, relate to 
definitely relativistic voltage or frequency ranges 
and are, therefore, not fully applicable. 

Values of J, and J,=I,, calculated by the 
methods described above are listed in Table II, 
along with the related values of V/Z? and v/v». 
Birge’s* 1941 values of the atomic constants have 
been used. Since the J’s vary approximately as 
the ninth power of e, precision here is not as 
superfluous as might be assumed. Indeed, an 
error in é as great as the discrepancy between 
today’s value and that of a decade ago could 
easily throw the calculated intensities off 
by much more than the errors of intensity 
measurement. 


*R. T. Birge, Rev. Mod. Phys. 13, 233 (1941). 


The quantity W, introduced in the fifth and 
last columns of Table II, is the total energy of 
continuous x-radiation in all directions. Inte- 
gration of J, of Eq. (2) over the 4m steradians of 
space surrounding the point of collision gives at 
once W=(8r/3)(J.+J,+J,], a function of » 
given in the units ergs per unit frequency interval 
per bombarding electron per atom-per-square- 
centimeter of target area. 

The numerical values in Table I] have been 
checked against each other by many graphical 
tests, and indeed, such checks were used at all 
stages of the calculations as a guard against 
errors of computation. Although each row of the 
table (except at v/y»=0 and some at v/vo=1) 
represents over 1200 elementary computing 
operations, it is believed that computing errors 
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Fic. 1. Screening constants @ (in cm) of the Sauter atomic potential distribution V,= (—eZ/r) exp (—r/a). 
Values shown by curve give best agreement with Fermi atomic potentials. 


and approximations have in no case impaired a 
calculated result by more than one percent. 


Low Frequency Limit 


The computations described above have been 
subject to the assumption that the bombarded 
nucleus is effectively unscreened. This does not 
preclude the application of the results to real 
nuclei,® fully clothed with screening electrons, as 
long as the acceleration of the bombarding 
electrons takes place inside of the surrounding 
electron screen. This is probably what happens in 
the production of most of the spectrum, but it 
definitely does not always happen when radiation 
of vanishingly small quantum energy is emitted. 
Such radiation requires the application of such 
slight forces to the bombardment electron that it 
may be produced by an electron passing a bare 
nucleus at a great distance. The collision cross 
section for radiation of zero frequency is very 
great, and one readily sees why the Sommerfeld 
theory predicts infinite intensity at A=. At 
this end of the spectrum, it is necessary to assume 
some screening if predictions applicable to real 
atoms are desired. 

We have, therefore, calculated J,,, at »=0 by 
the method of Sauter, who shows how to compute 
the radiation from an atom whose potential 

* A. L. Hughes, Phys. Rev. 55, 350 (1939), shows that 


atomic electrons play a negligible part in the elastic 
scattering of fast electrons by atoms. 


distribution is given by Vs = (—eZ/r) exp (—r/a), 
where r is a radial coordinate originating at the 
atomic nucleus. In the unscreened case (a=), 
Sauter’s method gives intensity values which are 
much too low in the high frequency region of the 
spectrum, but Elwert has shown that Sauter's 
method agrees with the Sommerfeld method in 


TABLE III. Radiation components at »=0 for 
screened nuclei. 














V Z I2X10" J, X10" V Z 12X10" J, X108 
28.43 4 .0584 .190 33.99 WwW 2.77 8.14 
44.40 5 0585 .199 48.94 36 2.78 8.45 
87.03 7 0586 217 83.40 47 2.78 8.93 

143.9 9 0586 230 131.4 59 2.78 9.34 
214.9 11 -0586 .240 201.2 73 2.79 9.78 
300.2 13 .0586 .248 299.1 89 . 2.79 10.2 
28.42 6 131 418 300.2 92 2.93 10.8 
50.54 8 132 451 
78.96 10 «132 476 33.36 35 3.80 11.1 
133.4 13 132 504 50.34 43 3.81 11.6 
202.1 16 132 527 82.36 55 3.81 12.2 
315.8 20 -132 544 133.4 70 3.81 12.8 
201.4 86 3.82 13.4 
31.58 8 .210 .673 230.5 92 3.82 13.6 
49.34 10 211 713 
83.40 13 211 758 50.04 47 4.57 2.99 
142.6 17 211 -803 
197.3 20 211 828 32.95 44 6.07 17.2 
308.4 25 211 864 49.64 54 6.08 18.0 
83.40 70 6.09 19.1 
50.04 13 351 1.17 144.0 92 6.10 20.3 
34.76 16 .763 2.39 33.36 46 6.56 18.6 
49.07 19 764 2.47 49.44 56 6.57 19.4 
84.93 25 .765 2.65 84.00 73 6.58 20.6 
139.1 32 -765 2.77 133.4 92 6.59 21.8 
206.7 39 -765 2.87 
300.2 47 766 2.97 33.15 58 10.5 28.7 
49.67 71 10.5 30.3 
300.2 70 1.70 6.36 83.40 92 10.5 32.3 
34.79 24 1.72 5.21 
50.80 29 1.72 5.46 
82.69 37 1.72 5.73 
133.4 47 1.72 5.98 
203.2 58 1.72 6.21 
296.0 70 1.72 6.44 
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Fic. 2. Radiation components from Table II at the long wave limit (v=0). (Note both axes are loga- 
rithmic.) J,>J, means negative polarization. 


its approach to infinite intensities at the low fre- 
quency limit. It seems reasonable therefore, to 
seek a proper value of a and apply the method 
at the long wave limit. 

The values of a adopted were those which 
brought the Sauter potential distribution into 
best agreement with the Fermi'® atomic potential 
Vr =Zeg/r over the region of r for which the 
Fermi expression is defined. Figure 1 shows the 
variation of a with Z. 

Equation (20) of Sauter’s" first paper may be 
specialized for y=0 and decomposed into 





e° f — 2h? 8m Vea? ) 
I,= log + 1 
cim( V/Z2)L8m Vea? h | 
- | 
I (7) 


 Am( V/Z?) 


1 h? 8mVeae ~ | 
‘(co era 
2 8m Vea? h? 


These values of J, and J, are not functions of 





” E. Fermi, Zeits. f. Physik 48, 73 (1928). The quantities 
eand ¢ are defined and tabulated in this paper. 

" There are two mistakes in the equation as printed: R® 
should be R?, and » in the denominator outside the braces 
should be 3. 


V/Z* alone but depend also upon V. This 
dependence is very slight for J,, a variation of V 
by a factor of 10 (V/Z*? remaining constant) 
affecting J, by much less than one percent. In the 
case of J,, the corresponding variation goes as 
high as 25 percent. The difference is classically 
understandable since the force between a nucleus 
and a remote passing electron is principally in the 
component lying along the impact parameter. In 
plotting spectral features as functions of V/Z* 
and v/v, we have in some cases to follow used 
mean values of the low frequency limit J, and J, 
within the narrow scatter described above. Such 
mean values are entered in Table II, and 
Table III contains a fuller display of all the 
separate values calculated from Eqs. (7). These 
easily calculated values (Fig. 2) are in every case 
readily acceptable extensions of our hard-won 
curves for finite frequencies. 


High Frequency Limit 


No short cut is available for computing the 
short wave limit radiation components. The 
elaborate methods described above are applicable, 
but several of the series must be carried out 
to tedious lengths. The second reference of 
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Fic. 3. Radiation components from Table II at the short wave limit (v=¥vo0). (Note both axes 
are logarithmic.) 


Weinstock indicates the special forms which the 
working equations take at this limit, but no over- 
all simplification results. We have, therefore, 
carried out the computations for only four cases 
at the limit. These have been supplemented by 
limit data obtained by extrapolation of results of 
related cases for which v/vo<1. Since the radia- 
tion components are functions of two variables 
(V/Z? and v/v), it was possible to arrive at the 
same result by two quite independent extrapola- 
tions. If two such extrapolations yielded results 
in agreement to within one percent, the result 
was accepted as correct. Figure 3 shows the 
variation of J, and J, at the high frequency limit. 
These curves are important in their relation to 
the mooted question of high frequency limit 
polarization. 


Approximate Representation 


Detailed inspection of Table II shows that the 
values of J, and J, are, to a rough approximation, 
inversely proportional to V/Z*. For J, the ap- 
proximation is close enough so that all the values 
within the range of our investigation may be 
stated with a mean error of 10 percent and a 
maximum error of 20 percent by the equation 


I,=3.0X10-"Z?/V. 


Near the short wave limit (v/v) =1), the variation 





of J, with V/Z? is not as simple, and the maxi- 
mum errors are found. Far from the short wave 
limit, approximate representations become quite 
satisfactory; all values of J, investigated at 
v/vo=0.2 are given by J,=2.6X10-"Z?/V with 
errors all under 2 percent. It should be realized 
that the word “error” as used here means only 
the discrepancy between the approximate ex- 
pression and the more rigorously calculated 
values of Table II; it has no relation either to the 
errors of the Sommerfeld theory or the errors of 
computation. 

The values of J, and J, are less easily summa- 
rized. They increase rapidly with increasing 
wave-length in a given spectrum whereas the 
general (but not invariable) tendency of J, is to 
decrease slowly. At v/vo=0.4 all values of 
I,(=I,) are given with errors under 6 percent by 
the expression 1.1 10-"'Z?/V, but elsewhere in 
the spectra, expressions of this type do not suffice. 

In view of the enormous amount of work in- 
volved in computing spectral characteristics from 
the Sommerfeld equations or any of their modif- 
cations, it seems worth while to consider a second 
approximate representation of their results. 
Using primes to distinguish these approximations, 
we find: 


I,’ X 108° V/Z?=0.252+a(v/v9—0.135) 
—b(v/vo—0.135)?, (8) 
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Fic. 4. Variation of I, (or J,) with position in the spectrum. Each curve pertains to a specific value of 
V/Z*. Reading downward from the top curve, these values are 0.00985, 0.01576, 0.0272, 0.0355, 0.0613, 
1axi- 0.1359, 0.493, 0.790, and 1.776 electrostatic units of potential difference. 
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) Fic. 5. Variation of J, with position in the spectrum. Each curve pertains to a specific value of V/Z?. 
it by Reading downward from the top curve these values are 0.00985, 0.01576, 0.0272, 0.0355, 0.0613, 0.1359, 
re in 0.493, 0.790, and 1.776 electrostatic units of potential difference. 


ffice. 

k in- where where 

from 

dif a=1.47B—0.507A —0.833, pw ota $1.21 y2—ys 
1.43y,—2.43 

cond b=1.70B—1.09A —0.627, | ellie 

sults. : j=(14+2h)y2—2(1+h)ys, 

=exp (—0.223V/Z*) —exp (—57V/Z?), 

—. Se anes ea k=(1+h)(ys+J), 


yi =0.220(1 —0.390 exp (—26.9V/Z?)), 
(9) y2= 0.067 +0.023/((V//Z2) +(0.75)), 
(v/v) +h ¥3= —0.00259+0.00776/((V/Z?)+0.116). 








B=exp (—0.0828 V/Z*) —exp (—84.9V//Z?). 











(8) Ii, .X10°V/Z?= —j+ 
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Fic. 6. Total intensity (ergs per steradian per unit frequency interval per bombarding electron per 
atom-per-cm*) emitted at 90° to the direction of bombardment as dependent upon position in the spec- 
trum. Each curve pertains to the value of V/Z? shown just above it. 


By these elementary expressions it is possible to 
compute the characteristics of continuous x-ray 
spectra with an expenditure of time and labor of 
the order of one percent of that required by the 
basic theory. Within the range covered by 
Table II, J,’, I,’, and J,’ differ from J,, J,, and I,, 
- respectively, by less than five percent at the 
worst and by about two percent on the average. 
The forms of Eqs. (8) and (9) have, naturally, no 
theoretical significance ; the equations are merely 
empirical fits of the data of Table II. 


Discussion of Table II 


Figure 4, a plot of data from Table II, shows 
the variation of the y (or z) component of radia- 
tion with its controlling variables. In any given 
spectrum the component J, (or J.) diminishes in 
magnitude as the short wave limit (v/vp=1) is 
approached. The diminution is most notable 
when targets of low atomic number are subject to 
bombardment by fast electrons. At the other end 
of the spectra, the Sommerfeld equations would 


have turned all the curves up toward infinite 
values. Introduction of screening in all calcula- 
tions for which v/v»=0 keeps the components 
finite but gives multiple values, since in this 
approximation the roles of V and Z do not always 
bring them together in the quotient V/Z*. This 
accounts for the frayed appearance of the left- 
hand ends of curves of Fig. 4. It is presumed that 
if it were practicable to take correct account of 
screening throughout the spectrum, the curves of 
Fig. 4 would, throughout their lengths, consist of 
multiple bundles converging closely with in- 
creasing v. } 

Figure 5 is a similar representation of J,. The 
insensitivity of J, to frequency (emphasized by 
logarithmic plotting) is the most conspicuous 
feature here. All points at v=0 have been calcu- 
lated by Eq. (7), but screening has little effect on 
the long wave limit values of J, and the splitting 
of the curves is too small to be shown. 

The total radiation emitted in the direction 
§=90° is given in accordance with Eq. (1) by 
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FiG. 7. Variation of total emitted energy W (in ergs per unit frequency interval per bombarding electron 
per atom-per-cm?) with position in the spectrum. The number immediately above each curve is the value 


of V/Z? to which the curve relates. 


Igo°=1.+T],. Since the direction is a favorable 
one for observations, we show in Fig. 6 the 
predicted total intensities. The simplicity of 
these curves is remarkable. The variations of J, 
and J, with frequency, being opposite in sense, 
have nearly neutralized each other and given for 
this direction of observation a set of almost flat 
and linear spectra. It should be realized that 
these ordinates are proportional to intensities 
within freguency intervals of equal width. Most 
methods of spectrum observation seek to measure 
a quantity proportional to the intensities in 
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Fic. 8. Polarization components of radiation originating 
at 0 by bombardment from the left and observed at B. 


wave-length bands of equal width,.and such 
spectra are not flat and horizontal in character, 
though they may be readily transformed into a 
form comparable with the curves of Fig. 6. 

Shifting the direction of observation to values 
of @ either greater or less than 90° tilts the curves 
of Fig. 6 so that their slight negative slopes be- 
come more accentuated and deepens their slight 
upward concavities. 

The total radiation energy W, from Table II, 
is shown in Fig. 7. 


POLARIZATION 


In Fig. 8 the origin of radiation is O, the origin 
of the coordinate system, and the direction of 
observation is along OB, a direction specified by 
the angle 6. The two vectors erected at B illus- 
trate the two components of radiation with 
mutually perpendicular directions of polarization. 
The component polarized with electric vector 
parallel to the y-axis of the coordinate system 
possesses intensity 7,; the component polarized 
in the perpendicular direction in general contains 
incoherent resolved contributions from both J, 
and J, and possesses the intensity J, sin’ @ 
+I, cos? 6. We designate the ratio of these com- 
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Fic. 9. Polarization of radiations emitted in a direction at 90° to the direction of bombardment. 
Numbers on the curves are values of V/Z*. Scherzer investigated the case V/Z?=0 and concluded that 
P 0° =0.6 for the entire spectrum. The curves above seem to be approaching this case as a reasonable limit. 


ponents by R=J,/(J, sin? 6+J, cos* @) and take 
as the measure of polarization the quantity 
P,=(1—R)/(1+R), an expression which is equiv- 
alent to Py=(1—D)/(1+D(2 csc? 6—1)), where 
D is Sommerfeld’s ‘‘depolarization ratio” J,/J,. 

From the data of Table II, the polarization of 
radiation emitted in the direction 6 =90° has been 
calculated for a number of values of v/vo and 
V/Z*. These calculations are presented in the 
curves of Fig. 9. The similarity of these curves 
shows that in the present theory, polarization is 
more a matter of position in the spectrum than of 
atomic number or bombardment energy. Only in 
the vicinity of the short wave limit does Pyo° re- 
spond markedly to variations of V/Z*. This 
response is in the sense that high V/Z? is associ- 
ated with high (though in no case complete) 
polarization. 

Near the low frequency end of the spectrum, 
radiation is produced at slight cost to the energy 
store of the bombardment electron. Close en- 
counters are not required for such emission, and 
the predominant force component acting upon 
the electron lies in the yz-plane rather than the 
x-axis. The classical consequence, negative polar- 
ization, is supported by the curves of Fig. 9. 


The curves of Fig. 9 present a somewhat 
tangled appearance, but they are in fact related 
in an entirely systematic manner, as a replot of 
the same data in Fig. 10 shows. The latter figure 
also emphasizes the relative insensitivity of 
polarization to V/Z? except in the region of low 
values of this quotient. 

To show the theoretical variation of polariza- 
tion with direction of emission, Fig. 11 is 
presented. Each curve here is a graph of Ps vs. 6 
for the particular value of v/v» appended to the 
curve and for the value V/Z? = 1.776 common to 
all five curves. Polarization of the forward 
directed radiation is always zero; as @ increases, 
the polarization may increase or decrease, de- 
pending upon v/v. Evidently for a given V/Z’, 
there is a value of v/vo for which the emitted 
radiation is substantially unpolarized for all 
directions of emission. The maximum (or mini- 
mum) of polarization at @=90° is very flat in all 
cases, and observers in this angular region need 
have little concern about the precision of their @ 
measurements. 

Alternative Methods of Calculation 


Sauter’s treatment of the problem made use of 
the first Born approximations. This method is 
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Fic. 10. Polarization of radiations emitted in a direction at 90° to the direction of bombardment. Numbers 
on the curves are v/vo. 


not appropriate unless Ze? /hv;, »<1, so the method 
is suited only to cases in which the bombardment 
electron moves with high speed both before and 
after collision. But such cases should properly be 
treated relativistically and with consideration of 
retardation of potential. 

Sommerfeld and Maue made an important con- 
tribution to the solution of the difficult problem 
of integrating the absolute squares of the matrix 
elements over all directions of electron departure 
when they showed that the integration necessary 
to the determination of the total radiation in all 
directions could be carried out exactly. The 
integral contains a hypergeometric function which 
they evaluated by an approximation which has 
the effect of giving spectral intensity values 
which are somewhat too low, as Elwert later 
showed. Elwert’s own method undoubtedly 
gives a better approximation, but it is not 
applicable at the short wave limit since it is re- 
stricted by the condition (Ze?/hv.) — (Ze? /hv1) <1, 
[(|m2| —|1|)<1 in Elwert’s notation] while at 


the long wave limit, it predicts infinite intensities 
exactly as does the Sauter treatment. 

Upon finding that the ‘otal spectral energy 
given by his own development exceeded that of 
Sommerfeld and Maue by a factor 


(a) 

| ny | T 

at the long wave limit and by almost the same 
factor at the other extreme, Elwert conjectured 
that the same factor might rectify the low in- 
tensity values of Sommerfeld and Maue in all 
cases. This supposition was strengthened by his 
investigation of its effect upon a single com- 
ponent of the matrix elements. 

We have tested this supposition by comparing 
certain of our own calculations (in all of which 
the hypergeometric function in question has been 
evaluated without error in excess of one percent) 


with corresponding intensities derived from the 
Sommerfeld-Maue approximation rectified by 
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Fic. 11. Polarization vs. 6 for V/Z*= 1.776 and for the values of v/v shown on the curves. Since in the 
theory here considered P¢ is a function of sin 6, the curves for the next quadrant will be symmetrical 


continuations of those shown here. 


Elwert’s correction factor. Elwert’s'® statement 
of the Sommerfeld result is equivalent to 








4nn, 

N,,?= { | 212+ 2?! log (1 — >) —4| mime! }, 
£072 
8rn, 

N,,?= { | my2—m?| log (1—£&) —4| myn! }. 
£o7m2 


The subscript s stands for Sommerfeld. Calcu- 
lation of J,, I,, by these equations and the 
Elwert correction factor is exceedingly simple. 
Figure 12 shows that the resulting values are in 
good agreement with corresponding values from 
our Table II for the small values of | 2! — || 
involved in these cases, In Fig. 13 a similar 
comparison is made for larger values of || 
‘—|m2|, and the agreement is much less good. 


Elwert’s proposal for rectifying spectral intensi- 








2G. Elwert, reference 3, Eqs. (34) and (35). 


ties obtained under the Sommerfeld approxima- 
tions is quite effective when his restricting 
assumption is satisfied. 


Efficiency of Continuous X-Ray Production 


From the radiation components of Table II, 
the efficiencies of continuous x-ray production in 
thin and thick targets may be computed. By 
thin target efficiency we mean the ratio of 
emitted x-ray energy (including all continuous 
spectrum frequencies and all directions of emis- 
sion) to the total energy loss suffered by bom- 
bardment electrons in passing through a sample 
of matter so thin that only a small fraction of the 
electron energy is dissipated. This efficiency is 
not simply a function of V/Z?, but of Z and V 
(or 8) separately. The method of calculation in a 
given case is to select from Table II values of W 
pertaining to a common V/Z? but to different 
values of v/vo and to perform a numerical or 
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Fic. 12. Radiation com nents for the case V/Z?=1.776. Solid curves from Table II. Broken 


curves by Elwert’s method. 
0.00829 < | m2| — |, | <0.167, and | n,| =0.160 


graphical integration yielding the total continu- 
ous X-ray energy per electron per atom-per-cm? 
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Fic. 13. Radiation components for the case V/Z? = 0.0604. 
Solid curves from Table H. Broken curves by Elwert's 
method, an approximation permissible only when |; |<C1 or 
|m2| — |m,|<C1. For these curves 0.112 < | m2! — |m| <0.545, 


and |»,| =0.867, 


agreement is found when | m2| — | 2, |<K1>>| ;|. For these curves 


in all frequencies from zero to the quantum limit. 
It then remains to divide this result by the total 
energy lost per electron per atom-per-cm’, a 
quantity depending upon 8 and Z. 

In seeking the most probable values of this 
quantity, it seems best to adopt a combination of 
theoretical and experimental results. Williams" 
has shown that the observations of White" and 
Millington upen the stopping of fast 8-particles 
by mica agree with the classical theory of Bohr'® 
with respect to the functional variation of stop- 
ping power with 8, though not as to the absolute 
values of stopping power. Williams recommends a 
practical formula consonant with these facts 
which may be put in the form :'* 

8 E. J. Williams, Proc. Roy. Soc. A130, 310 (1931). 

4 P, White and G. Millington, Proc. Roy. Soc. A120, 701 
“ N. Bohr, Phil. Mag. 25, 10 (1913). 

‘6 An error in the statement of the energy-loss formula in 
Williams’ Eq. (18) has been corrected above. We have also 
incorporated a suggestion of D. L. Webster that the right- 


hand side of Williams’ formula be multiplied by 22/A 
where A is the atomic weight, 
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Fic. 14. Rates of loss of kinetic energy of electrons of 
speed 8 in the four substances indicated by atomic num- 
bers near the curves. Ordinates are in ergs per atom-per- 
cm? X 10” of stopping material. These curves are adapted 
from E. J. Williams’ discussion of the theory of Bohr. 


9 


Energy loss in ergs per electron per atom-per-cm? 
| =5.62Z8-"-*X10-**. (10) 


This formula takes no account of the losses by 
those few collisions in which an electron loses 
more than 1500 kev; its applicability above 
8=0.52 has not been demonstrated nor has it 
been tested with elements of high Z. Nevertheless, 
we adopt it for all speeds above 6 =0.5 and for all 
elements here considered. 

For low values of 8 (up to 8=0.30) we draw 
upon Williams again. By cloud chamber observa- 
tions of x-ray photoelectrons in oxygen, hydro- 
gen, and argon, he established that in this energy 
region the Bohr theory predicts rates of energy 
loss about 60 percent too high, but correct in 
respect to their functional dependence upon Z 
and 8. Assuming that the same ratio holds for all 
elements we multiply Bohr’s expression by 0.6 
and obtain: 


KIRKPATRICK AND L. WIEDMANN 





Energy loss in ergs per electron per atom-per-cm? 
=4.07PZp X10. (11) 


In this equation, P is a dimensionless function of 
Z and V involving the geometric mean of the 
natural frequencies of all electrons in the atom 
of the stopping material. P was calculated for 
Z =13, 47, 70, and 92 by use of frequencies from 
the x-ray v/R term values found in Siegbahn’s 
Spektroskopie der Réntgenstrahlen (1931). Results 
of these calculations appear in Table IV. 

To obtain electron energy losses in the region 
of energies intermediate between 20 and 80 kev, 
we have simply filled in by graphical interpolation 
in the manner illustrated by Fig. 14. It is 
gratifying that the calculations in the separated 
ranges admit of a smooth connecting curve seg- 
ment for each of the four elements used in these 
efficiency calculations. We do not take this to 
mean that these curves are entirely correct ; they 
may be wrong by many percent, but they would 
seem to be as good approximations as the litera- 
ture of the subject affords at this date. Our 
efficiency conclusions will be subject to immedi- 
ate amendment upon the appearance of more 
reliable energy-loss information. Figure 14, show- 
ing the energy-loss rates which we have adopted, 
may be used as a basis for any such adjustment. 


TABLE IV. Electron-loss rates from Eq. (11). 











Loss 
Kilo- rate 

Zz 8 volts P x10 
13 0.10 2.57 7.47 39.5 
13 0.15 5.84 8.38 19.7 
13 0.20 10.53 9.97 13.2 

13 0.25 16.75 10.77 9.12 

13 0.30 24.66 11.43 6.72 
47 0.10 2.57 5.66 108.2 
47 0.15 5.84 7.12 61.0 
47 0.20 10.53 8.15 39.0 
47 0.25 16.75 8.96 27.4 
47 0.30 24.66 9.61 20.4 
70 0.10 2.57 5.46 155.5 
70 0.15 5.84 6.92 87.6 
70 0.20 10.53 7.95 56.7 
70 0.25 16.75 8.76 39.9 
70 0.30 24.66 9.41 29.8 
92 0.10 2.57 4.91 183.8 
92 0.15 5.84 6.42 106.8 
92 0.20 10.53 7.40 69.3 
92 0.25 16.75 8.21 49.2 
92 0.30 24.66 8.86 36.9 
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Fic. 15. Efficiency of production of continuous x-rays by thick targets of atomic numbers shown 
beside the curves when bombarded by electrons at the stated kv values. Curves are based on 
electron energy-loss rates of Fig. 14 and theoretical continuous spectrum energies derived from 
Table Il. Approximate proportionality of efficiency to bombardment energy is shown. 


The thin target efficiencies are given in Table V 
along with thick-target efficiencies now to be 
explained. 

If E represents the varying kinetic energy of a 
typical electron slowing down from Eo to zero in 
a thick target, then R, the output of x-ray energy, 
is given by 


Eo 


(dR/dE)dE, 


TasBLe V. Calculated efficiencies of continuous x-ray 
production in thin and thick targets of atomic number Z. 
For thin targets kev denotes the energy of bombardment 
at the point of production. For thick targets it denotes the 
energy at incidence upon the target. 














kev Z=13 Z=47 Z=70 Z=92 

15 5.3 23.0 

25 7.9 34.0 50.0 69.0|\thin target effi- 
40 13.0 53.0 80.0 110.0/ cienciesXx10* 
90 27.0 110.0 160.0 210.0 

16.8 34 14.0 ) 

24.7 44 18.0 280 37.0) 4-4 in 
46.5 7.6 32.0 48.0 64. ceacian >< 10" 
79.0 13.0 51.0 79.0 100.0) “ences 

90.0 15.0 58.0 880 110.0 








where the derivative in parentheses is just the 
thin target efficiency. It is convenient to perform 
the graphical integration with respect to B,"in 
which case we have 


80 
R= f (dR/dE)(dE/dp)dg. 
0 


The second of the derivatives in the integrand is 
obtained by differentiating the relativistic ex- 
pression for electron kinetic energy in terms of 8. 
The integration has been carried out for twenty 
representative cases, and the over-all thick target 
efficiencies of continuous x-ray production have 
been obtained by dividing R by the initial 
electron energy. These efficiencies are shown in 
Table V, and Figs. 15 and 16. 

From direct measurement of the x-ray and 
cathode ray energies Beatty'’ found thick target 
efficiencies to be proportional to the product of 
target atomic number and bombardment voltage. 
The many efficiency measurements made subse- 
quently to Beatty’s time were summarized by 


7 R. T. Beatty, Proc. Roy. Soc. 89, 314 (1913). 
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Fic. 16. Efficiency of production of continuous x-rays by thick targets of atomic numbers shown by the 
ordinate scale when bombarded by electrons at the kv values stated at the ends of the curves. Curves are 
based on electron energy-loss rates of Fig. 14 and theoretical continuous spectrum energies derived from 
Table I]. Approximate proportionality of efficiency to atomic number is shown. 
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Fic. 17, Calculated efficiency of production of continuous 
x-rays in thin targets (upper curve) and thick targets 
(middle curve) of atomic number Z bombarded by electrons 
accelerated by potential kv. Approximate linearity of thick 
target curve supports Beatty's rule. Failure of points to lie 
precisely on the curves shows that efficiencies as here 
calculated are not strictly functions of the product Z kv. 
The lower curve is a consensus of the experimental results 
of several investigators of thick-target efficiency. 


Compton and Allison'’ in the equation Efficiency 
=1.110~-*Z kv, asummary which these authors 
regarded as ‘probably correct within about 20 
percent.” 

Figure 17 shows that the efficiencies of Table V"* 
agree approximately with Beatty's rule and that 
they lie about 20 percent above the values of the 
consensus of Compton and Allison. 

For any product Z kv the calculated thin 
target efficiency is about twice as great as the 
thick target efficiency. The thin target points in 
Fig. 17 are tolerably represented by the straight 
line 

Thin target efficiency = 2.8 X 10~®Z kv. 


Calculated efficiencies for both thin and thick 


18 A. H. Compton and S. K. Allison, X-Rays in Theory and 
Experiment (D. Van Nostrand Company, Inc., New York, 
1935), p. 90. 

1 The efficiencies of Table V lie somewhat below the 
values given earlier by one of the present authors (P. 
Kirkpatrick, Phys. Rev. 66, 156 and 161 (1944)). The 
original calculations were based upon the energy-loss rates 
of Eq. (10) for low voltages as well as high. The values of 
W at v=0 in the previous work were obtained by graphical 
extrapolation instead of direct calculation. Changing the 
method of calculation in these two respects produced 
efficiencies more precisely proportional to Z and kv and, on 
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targets tend to fall below linearity at high values Dr. Robert Weinstock assisted the authors 
of Z kv, but the reliability of the calculations is through numerous discussions in the early stages 
in question here since the potentials involved are of the work, when his interpretation of his own 
outside the proper range of non-relativistic papers contributed most helpfully to the de- 
theory. velopment of the computational procedure. 
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The Relativistic Correction in the Meson Theory of Nuclear Force 


Nine Hu 
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In the present paper a systematic investigation is made on the relativistic corrections of the 
usual non-relativistic meson theories of nuclear force by expansion with respect to the dimen- 
sionless operator V/M. The mixed theory initiated by Moller and Rosenfeld is considered 
since it is the only satisfactory theory in the non-relativistic region. It is found that the inad- 
missible singularity which has been removed in the mixed theory in the non-relativistic region 
reappears in the higher approximations. This means that the extension of the mixed theory 
to the relativistic region as done by Moller and Rosenfeld in explaining the quadrupole moment 
is not justified. 


I. INTRODUCTION 


HE well known fact that the deuteron nucleus possesses an electric quadrupole moment is an 
indication of the existence of the tensor force in the nuclear interaction. The expression for 
this force derived from a single type of meson field, however, contains the inadmissible singularity 
of the type 1/r*, which has to be cut off at an arbitrary radius. It has been shown by Moller and 
Rosenfeld! that if the nuclear interaction is assumed as due to a mixed field of pseudoscalar and 
vector mesons of the same mass and suitably chosen coupling constants, then the tensor force 
disappears with its inadmissible singularity in the non-relativistic region. This provides a solution 
of the wave equation without resorting to the cut-off procedure. According to this theory, the quad- 
rupole moment can only be a relativistic effect owing to the vanishing of the tensor force in the 
non-relativistic region. Mgller and Rosenfeld have actually derived an expression for the quadrupole 
moment by considering the relativistic interaction of the first order. Schwinger’ later suggested that 
if the mass of the vector meson is assumed to be larger than that of the pseudoscalar meson, then 
the tensor force reappears with an admissible singularity of the type 1/7. This assumption is also 
in agreement with the hypothesis that the vector meson is highly unstable and responsible for the 
8-disintegration of the nucleus. 

Recently Jauch and the present author* calculated numerically the quadrupole moment according 
to Schwinger’s assumption, and found that the result is much smaller than the experimental value 
even for a very large mass ratio for these two types of mesons. This conclusion is quite independent 
of the coupling constants obtained from the scattering problems. Parallel to this calculation, Hulthén‘ 
also evaluated the quadrupole moment using the formula obtained by Moller and Rosenfeld. The 


'C, M@ller and L. Rosenfeld, Kgl. Danske Vid. Sels. Math.-Fys. Medd 17 (1940). 
? J. Schwinger, Phys. Rev. 61, 287A (1942). 

3]. M. Jauch and N. Hu, Phys. Rev. 65, 289 (1944). 

‘LL. Hulthén, Arkiv for Mat. Astr. Och. Fys. 29A, No. 33 (1943). 
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result obtained by him is in reasonable agreement (judging from the error admissible by the approxi- 
mation) with the experimental value. This would justify Mgller and Rosenfeld’s original view that 
the quadrupole moment is essentially a relativistic effect and therefore their explanation would be 
nearer to the fact than Schwinger’s explanation. 

The original derivation of the formula for the quadruple moment by Mller and Rosenfeld is, 
however, not altogether consistent. In its derivation, the relativistic force is added as a perturbation 
term to the unperturbed static wave equation without considering the other terms of the same order 
of magnitude arising from the non-relativistic approximation. In the present investigation we 
retain these terms consistently and find that this non-static force gives rise to a tensor force also of 
the inadmissible type 1/r*, which in turn requires further adjustment of the coupling constants to 


make it vanish. Consequently, the relativistic effect, if made regular, cannot make any contribution 


to the quadruple moment. 

It should be noted that the inadmissible singularity of the type 1/r? also exists in the equations 
used by Mgller and Rosenfeld. No serious attention, however, has been paid to this since the first 
perturbation leads to Convergent and unique results. Application of the perturbation method in our 
case also gives convergent and unique results in the first perturbation. These results can only be 
considered as provisional and further investigation is needed to find out their significance, as we 
know that the equations give no stationary states when solved rigorously unless additional assump- 
tions equivalent to a cut-off are introduced. 


Il. THE NON-RELATIVISTIC APPROXIMATION ON THE DEUTERON WAVE EQUATION 


The rigorous relativistic interaction potential obtained by Kemmer,' Bhabha,*® and Moller- 
Rosenfeld! for a mixture of pseudoscalar and vector mesons is given by 


V=F+W-+K, 


1 
F= (1020 gst +g(o'e®) ~<a (oor ny (aor boi, 
ue 


W= (e070) (pvrgin — pe ) eT gellig tt x. gattngth 
K a 
+ (61a + pps) (@ Ko) ]]} VOR), 
(1) 
K= (x4) — fips ps pi py —g 1291 pr (@ Me?) 


9 


£1 £2 ae ‘. 
=p p (@ PV) (eV) ——— pi pr ps ps (@ YT) (@V) a(r), 
xy od 


—pr 


VO=—-Y2=V, 





7” 4nr 
where F gives the static interaction. W and K are the non-static interactions linear and bilinear 
in pi” and p2” (v=1, 2), respectively. u is the mass of the mesons (we are now considering the 
case that the pseudoscalar and the vector mesons have the same mass), 1, p2, p3 and @ are the usual 
Dirac matrices. (v)=1, 2 refers to the first and the second nucleus respectively, and r“? and r® are 
the isotopic spin matrices. f; and f: are, respectively, the spin-independent and the spin-dependent 


5 N. Kemmer, Proc. Roy. Soc. 166, 127 (1938). 
*H. J. Bhabha, Proc. Roy. Soc. 166, 501 (1938). 
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coupling constants for the pseudoscalar meson, and g; and ge are the same constants’ for the vector 
meson. The unit c=f=1 is used in the present paper. 
The relativistic wave equation is then 


1a 
wr ae Mx eV) + py (Pe (2720) M(p3°+ p3)+V v(x, x) = 0, (2) 
L 1 


where V is given by (1). ¥ is the spinor with sixteen components formed from the direct product 
of the two four-component spinor representations of the two nucleons. We may conveniently 
represent y by Was (a, 8=1, 2, 3, 4), where the first index goes only with the operators referred to 
the nucleon (1) and the second index goes only with those referred-to the nucleon (2). 

It is seen from (1) that the tensor force of the form (eV) (e¢@V®)a(r) gives rise to an inad- 
missible singularity of the form 1/r* which gives no stationary solution for the deuteron wave 
equation. If we put g:°=g,.*, then we get the well-known Mller-Rosenfeld mixture of which the 
tensor force with its inadmissible singularity disappears entirely from the potential F in the non- 
relativistic region. It should be remembered that Mgller and Rosenfeld’s theery was originally 
proposed for the static region only. In the relativistic region, as we can easily see from (1), the inad- 
missible singularity of the type 1/r* still remains. Therefore it is meaningless to solve the relativistic 
equation (2) rigorously for the Mgller-Rosenfeld mixture, as the latter gives no better service than 
the other simple theories. The situation will be different, however, if we try to solve (2) by the 
non-relativistic approximation of which the first approximation is the static case. Then it is quite 
justified to start with the mixed theory of Mgller and Rosenfeld which alone gives regular singularity 
in the first approximation. The second approximation would correspond to the region in which the 
M@gller-Rosenfeld expression for the quadrupole moment was derived. 

To carry out the approximation, we follow the well known procedure’ of expanding y in powers 
of the dimensionless operator V/M where 1/M is the Compton wave-length of the nucleon. We put 


(Wir, Wiz, Wai, H22) = (1, G2, Gs, gale, 
(Ws1, Ws2, War, a2) =(x1™, x2, x8, xa Emme, 
(Wis, Yrs, Vos, Hoa) = (x1, x2, x8, xa Orme, 
(Was, Waa, Was, Vas) = (G1, 2, bs, Gale™. 


Substituting (3) into (2), we obtain 


(3) 


10 
-=t4- oP) 4 “(0V) x4 Voge Voxnx® + Voginx® + Vey =0, (4a) 
1 


1 dx 
merece ding MeV) y4- (o¥)9— 2Mx 
t 1 
+ Vym—eet Vmax + Vy aoyiox © + Vy «gd =0, (4b) 


Ae 1 
"Tal sal dal? (00) 9— 2Mx 





4 
+ Vycnget Vyinyox + Vy oyinx + Vy igh =(), (4c) 


1 d¢ . 
-- 4. te Y)y@ yw “(6 29 2))¥ 0) —4MO+ Vepot Voxwox™ + Vexenx™ + Vege =0, (4d) 
1 1 


’The nomenclature that g, is the spin-independent coupling constant is only justified in the non-relativistic region. 
In the relativistic region, as we can see from (1), the spin-dependent interaction terms also contain g:. 
8 W. Pauli, Handbuch der Physik (Verlagsbuchhandlung, Julius Springer, Berlin, 1933), Vol. 24, Pt. 1, p. 237. 
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where ¢, x“, x, @ are the four-component wave functions defined by (3). o and o® are 


respectively the direct products «XJ and I Xe, @ being the usual Pauli matrices and J =(51): 


V,«¢ is the submatrix of V that transforms a ¢ into a x“ and Vy i, Vyaxcm, +++ have similar 
meanings. We can easily verify from (1) that 


Voe= x(x) FV y(y (i = Veo =F, 


Vox = xHe = x(neo = Vex» =Gi, 


(S) 
Vexi = x(vde = xg = Vexa) =A, 
Voo = Vee = xxi = x(oxa) =Ky, 
where 
1 
Fim(se) ett ertoe®) —e — £2 aa my(oovoy(oove lo(), 
7 
Gy=(r\ 2) AA cg E900 440% X0)]] -VOa(r), 
7 
(6) 
HMy=(1r%r {Ae a —- e Ma -+ei(o®X0%)]] -V5(0), 
7 


2 


gi" 
K,=(1 7) — fites Pes —g1(0 a) += (69) (0 Y) +e 5(eg® (0) (2) a(r), 
ld lad 


where f.?=g,.” since we are now considering the M@ller-Rosenfeld mixture. The quantities ¢,”? are 
defined in Table I, where the first column gives the value of e,”’ when Fi, Gi, Hi, or Ky is 
applied to yg, and similar meanings for the the other columns. 

Before introducing the non-relativistic approximation, we shall first estimate the order of 
magnitude of the different terms in (4) with respect to the dimensionless operator V/M. The expec- 
tation value of the operator V is roughly of the order 1/ro, where ro is the radius of the deuteron 
nucleus. It follows that V/u is of the order unity since n.~7r,~'. On the other hand V/M is a quantity 
much smaller than unity (roughly of the order n/M~0.1). Therefore V/M is a suitable parameter 
for a perturbation approximation. Furthermore, in the derivation of (1) a perturbation expansion 
with respect to the parameters g./(41)! and f./(4)* (a=1, 2) has already been used. The order of 
magnitude of g,?/4r and f.?/4# determined from the deuteron problem is also found to be very 
small (roughly ~0.1). Therefore in the following approximation we shall also count g.*/4m and 
fa*/4m as small quantities of the first order. In other words, whenever terms of the mth order in 
V/M are neglected in the part proportional to g,’, f.?, or faZa, only terms of (n+1)th order are 
neglected in the part independent of the interaction. In the following we shall refer the order of 
magnitude only to the interaction part of our equations. The present method of expanding y’s in 
powers of V/M is essentially the non-relativistic approximation since V/M applied to ¢ gives the 
velocity of the nucleon measured in the unit c=1. We shall see later that the usual non-relativistic 


wave equation can only be obtained when terms 











of the order V/M are neglected in (4). TABLE I. 

Since only (4a) of (4) does not contain the " = aaa 
term of the form M we may consider ¢ as zeroth 3) 1 ct 1 aa 
order of magnitude. We can easily verify from €,) 1 1 -1 -1 
(4b, c, d) with the help of (5) and (6) that x i : “ a = 








and x®) are of the order V/M, and ¢ is of the 
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order V?/M?. Satta terms of the order V?/M? in (4a) and of the order V*/M? in (4b, c, d), we 
obtain 


1d¢ 
aan oats Hemym)y a4! slo) x + Vewet Vexox™ + V,xox® =0, (7a) 
¢ 

1 
2Mx =—(e9V) 94+ Vy weg, (7b) 

1 

1 
2Mx® =—(eV ®) 9+ Vy cre, (7c) 

4 
4Mo= Vege: (7d) 


Substituting (7b) and (7c) into (7a), we obtain 


1d¢ 


1 
—- 2 yw $V" ]0+ Vero t——{(0V)(Viurge) + Verne OV)e 


+(¢2V®)(Vy rey) + Vexn(oV) 9} =0, (8) 
where from (6) 
Voo= (rrp +2:%(0e) Jar), 





—_ =(r40f Fafa cay 4 Bubs io Eom xo)| -VM5(r), 
KB ke os 
Vycrg =(17 7 | gt — £183 1) FE (e xe]. Var), 
rm Me rm (9) 
Vx (rr){ fife. 209 4 5, “i 2) — 2? 902) x g(t) -V25(r), 
bb a be 
Vytng = (477 ogc — 6182 (2) 8188 9 (2) g(t) V(r). 
rn M L 
If we further neglect the terms of the order V/M in (8), then we have 
ao oe vay’ +9 (2) . + V, =0, 10 
i at 2M ¢ ee? = (10) 


which is just the non-relativistic wave equation for the deuteron with the static interaction potential 
Vo» of Mgller and Rosenfeld. It should be noted that [(1/M*)V*a(r)] is no longer a small quantity 
of the order m when r is very small, since @(r) behaves as 1/r for small r. Consequently the neglect 
of the terms of the form [(1/M*)V*a(r) ]g in passing from (4) to (8) and (10) is justified only when 
¢ is also very small for small r, so that the contribution of these terms is still negligible. This fact 
can be utilized as a test of the validity of our approximation, since ¢ can be determined completely 
from the wave equations. Physically this means that the two nucleons must always remain suf- 
ficiently apart from each other in order that the present non-relativistic approximation can be used. 

Equation (8) is the approximate wave equation for the deuteron accurate up to the first order 
in V/M. Substituting (9) into (8) we obtain after some calculation 


1d¢ 
-52- Save" +702) “et V. eet Pe =0, (11) 
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with 
1 
Pe= Fy Or aha erg) OPV (COV) a]-9 
& 
1 
tele este (eM +0 VOGXVY 9) +(e +e), V2OaXV¢) ] 
& 
1 
pre este. 1 + (oM0™) Inte o, (12) 
bb 


where we have replaced V’a by u’o in the third term of (12), since V?@—y?a =4(r), the delta function 
in the interaction being neglected as usual. It can easily be recognized that the singularity of the 
first term of (12) is of the inadmissible type 1/r*. To investigate the nature of the singularity of 
the second term, we put it in a slightly different form; using the center-of-mass coordinate system 
of the deuteron in which V g= —V®g, 


1 
Mi eedle te™, VOSXIV og) + (e+e, V2MOaXV¢) ] 
M 


I (o+0%, ‘xve) Pe ee. co (e+e), Ly), (13) 
Mu dr r Mu r dr 


where L is the angular momentum operator defined by 
L,=y(0/dz)—2(0/dy), L,y=2(0/dx)—x(d/dz), L.=x(d/dy) —y(d/dx). (14) 


It has been shown by Kemmer’? that the form of ¢ for the ground (triplet) state of the deuteron is 











g=Agiot Beginz, (15) 
_ (bu a AiO _( [m(m+1)}P¥or, — —[(m+1)(1—m)} Yo) 1 Ro(r) i 
ee Nd oD NE +¢-m)(1—m) PY", Cm(m—1) PY v2 or 
a Yiu ea) wo(r) _ [(2—m)(3—m) }'Y2"""', [(2-+m)(2—m) }}Y." 1 Riv) (17) 
Fe Ne v2) 2? oe AE (24+m)(2—m) PY",  [(2+m)(3+m) PY) 25 7’ 


where A and B are two arbitrary constants satisfying A?+B*=1, / is the orbital quantum number’ 
and m=0, +1 are the quantum numbers for the component of the spin in the z-direction. Y," is the 
normalized spherical harmonics, ¥)o”" and ¥).” stand for the matrices in (16) and (17) respectively, 
and R,; and R: are two normalized radial functions. Substituting (15)—(17) into (13) and noticing 
hat Lyg,.»=0, we find after some reduction 


1 
Wer gesllo +e, VS XVM g) + (eM +e, VSXVy)] 
e . 


2 ( (1) (2)) l de B 12 
= — T ———- 9 = 
Ma wnt bee 


1 
(11) guge( + Jew Boras (18) 


r? 





Therefore the second term of the right-hand side of (12) also gives rise to a singularity of the inad- 


*N. Kemmer, Helv. Phys. Acta 10, 47 (1937). 
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missible type 1/r*. After some reduction on the first term, (12) becomes 





3 
Pperm—tety ©) fufeteee)[on)(on) —H(o%e)}(—4~ - te ewe 
r 


4nh 


' 4? 
4+-— — (1) 7’) (1) nl ap. djs 
— Vfifet+egige)- (ee - e ete 





(7978 avee( — >= “)e “ Boins 


—yr 


(772) gf 1+(e Pe) u2—-~ (n=r/r). (19) 
r 





 & Mu 


Substituting (15), (17) into (11) and (19), using the relations (o@e®) =1, and (r“r@) = —3 for 
the ground state of the deuteron, and equating the coefficients of 9)" and of J).", respectively, in the 
result, we obtain finally 


1 d*u eo 
| eat gs" 








ae 


Mu 


Tr 


(fifet eet <—u20, (20) 
ae emmnoomes ge——-u =0, (2 
4xM eed eT al 








1 {dv 6v 
-—( fl [2 Eo+3° an en b— 22 


dr? r* mr 


(fife+ (5 +24") 
9 ” —_ — — “are 
4nMu vial Pr of re 


—yr 








iia 
Sa iht ees) (5 es ee ry 


(++) 0, (21) 
om « 2 oumes a= e~*ry+ 9 =U, 
ii ry? 4r J ee 


f (u?+-v?)dr = 
0 


(20) and (21) are the two radial equations for u and v derived from (11). It is seen that the inad- 
missible singularities appear with coefficients proportional to either (gigo+/if2) or gige. To remove 
these singularities the only way is to put gige=fife=0. This is, however, equivalent to putting 
fi=gi=0, since fz and gs have already been fixed in Moller and Rosenfeld’s original non-relativistic 
theory. 

With the reappearance of the inadmissible singularity in the relativistic approximation, our 
wave equation (11) will give no stationary states in the rigorous solution unless some assumption 
equivalent to a cut-off is introduced. This difficulty is the same as that encountered by the single 
type of meson theories before the introduction of the M@ller-Rosenfeld mixture. A cut-off is thus 
inevitable, and the application of the mixed theory is no longer justified as it gives no better service 
than the other simple theories. On the other hand, one can remove the inadmissible singularities 
by putting f:=g:=0, but then the tensor force disappears entirely from (20) and (21) and the 
observed quadrupole moment still remains unexplained as in the non-relativistic theory. It seems 
that the situation will remain the same in all the higher non-relativistic approximations, namely, 
whenever we obtain a tensor force it is always inherent with the inadmissible singularities. Further- 
more, the removal of these singularities will then be impossible since all the interaction constants 


have already been fixed. 


bl Mu 











where «=AR)(r), v= BR,(r) so that 
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Without considering seriously the difficulties mentioned above, Moller and Rosenfeld, and later 
Moller himself,!® tried to solve their equations corresponding to (11) which also contain the inad- 
missible singularities, by the perturbation method, considering the relativistic terms as the per- 
turbed terms. They obtained convergent and unique results for the quadrupole moment in the first 
approximation. Their results thus obtained can only be considered as provisional. Some troubles 
are therefore bound to be encountered in the higher approximations in such a way that the per- 
turbation cannot converge. To avoid this divergence, we have again to introduce the cut-off pro- 
cedure as is inevitable in the rigorous solution. However, it seems quite likely that the arbitrary 
radius of cut-off will only influence sensitively the higher approximations where the results diverge 
Originally before the introduction of the cut-off. This would mean that the results obtained by 
Mller, Rosenfeld, and Hulthén still give the principal part of the quadrupole moment, which is 
independent of the arbitrary radius of cut-off. Further discussion of this point will be the topic of 


another communication. 
For the sake of completeness, the same perturbation method as used by Moller and Rosenfeld 


has also been applied to solve (11). We find 
Et Eo 
(EA, l=2|P| Eo, 1=0) = f ei-sPoinads 


2v2 "75 3u pe y ; 
= ———(fifot+g1g2) (++ +- Jere Rote) ReM(rdr, (22) 
4nMu o \P Pr fF 


where Eo is the deuteron binding energy and E is the continuous positive energy. The quadrupole 
moment is given by (cf. Mgller and Rosenfeld’s paper) 


ats ll pare deine didn st (23) 


5 Jo Eo—E 





v2 ¢* (l=0, Eo|r?| E, 1 =2)(E, 1 =2| P| Eo, 1=0)+conj. . 





Substituting (22) into (23), we find that the result is identical with final expression (123) for quad- 
rupole moment in Mgller and Rosenfeld’s paper. Therefore the quadrupole moment obtained from 
the present theory by the perturbation method is the same as that evaluated by Hulthén. 

In conclusion, the author wishes to express his sincere thanks to Professor W. Pauli for suggesting 
the present problem and for his invaluable help during the investigation. 


°C, Moller, Kgl. Danske Vid. Sels. Math.-Fys. Medd 18, 6 (1941). 
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The Azimuthal Variations of Cosmic Radiation at 22° Latitude 
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Four triple-coincidence cosmic-ray telescopes directed at a common zenith angle of 60° were 
mounted on a turntable. Each telescope represented a vertical plane 90° from that of its neigh- 
bor. In each was inserted” 10.2 cm of lead. By 180° reversals each telescope was interchanged 
in its position with the one directly opposite. Countings were made at settings of the table for 
every 10 degree of azimuth angle. Azimuthal variations in the North-West and North-East 
quadrants check qualitatively those theoretically predicted by Hutner. The data are still not 
sufficient to warrant quantitative comparison. It has been found that the East-West asymmetry 
does not coincide with magnetic east-west plane, but occurs at 280°-100° plane. 


R. HUTNER!' made a quantitative study 

of the azimuthal effect of cosmic radiation 
(variation of cosmic-ray intensity at a given point 
on the earth’s surface and fixed zenith angle but 
for different azimuths), using a constant zenith 
angle of 60 degrees, at a geomagnetic latitude of 
20° north. This geomagnetic latitude was chosen 
because the penumbra? (the reader is referred to 
Shremp’s paper for reference for details about the 
penumbra) which makes possible the determi- 
nation of the energy spectrum of the primary 
cosmic rays makes a significant contribution at 
intermediate latitudes (between 10° and 40° 
geomagnetic latitude). The author undertook to 
perform an experiment at Lahore (geomagnetic 
latitude 22° N) to find the predicted azimuthal 
effect for a fixed zenith angle of 60°. A prelimi- 
nary report® on the experiment appeared in the 
Physical Review in 1941. This is therefore a 
progress report on the experiment. 

In the present experiment four cosmic-ray 
telescopes were used, each consisting of three 
triple coincidence G-M counter tubes with lead 
filter of 10.2 cm. The distance between the 
centers of the first and the third tubes was 36 cm 
in each case. The solid angle range covered by 
each telescope was 8.07° in the vertical and 
38.87° in the lateral planes. In the previous ex- 
periment® the solid angle range covered by each 
telescope was 14.5° in the vertical and 65° in the 
lateral planes. The telescopes were mounted on 
the four edges of a square table, which was 


'R. A. Hutner, Phys. Rev. 55, 614 (1939). 
?E. J. Shremp, Phys. Rev. 54, 157 (1938). 
+P. S. Gill, Phys. Rev. 60, 153 (1941). 
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rotated back and forth through an angle of 180° 
by hand at frequent intervals. The frequent 
interchanging of the opposite telescopes made it 
possible to compare the cosmic-ray intensities 
from opposite directions. The azimuths were 
changed by steps of 10° each time. Each of the 
four telescopes swept the whole sky for the fixed 
,zenith angle at intervals of 10° at least twice, and 
some azimuths were covered more than twice. 
The two or more readings of one telescope at the 
same azimuth, though taken about a month and 
a half apart, do not differ from each other by 
more than their statistical fluctuations and thus 
showed the systematic working of the telescope 
for the whole experimental period. These readings 
were taken between October 31, 1941, and 
February 21, 1942. 

Compared with the solid angle range of the 
earlier experiment, that of the new experiment is 
reduced giving therefore a resolving power three 
times as great. The number of counts per unit 
time for each new telescope should therefore be 
one third of the number of counts per unit time 
in the earlier experiment. The only other differ- 
ence between the new telescopes and the old was 
that, among the four, some of the individual G-M 
tubes had been interchanged. Another rough 
check on the systematic working of the tubes was 
provided by the fact that the ratio of the counts 
per day in the new experiment to that in the old 
was 1:3.12. 

The apparatus was installed in a special labo- 
ratory constructed for the purpose on top of the 
new physics laboratory of Forman Christian 
College. The roof of the observatory was made 
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of a single sheet of galvanized iron and was thus 
virtually open to the sky. 


INTERPRETATION OF DATA 


Curves A, B, C, and D of Fig. 1 represent the 
number of triple coincidences per hour obtained 
by telescopes A, B, C, and D for each 10° interval 
of azimuths for the fixed zenith angle of 60°. 
Each telescope made a complete turn of ?the 
horizon, except for the azimuths where one of the 
amplifying tubes or one of the power tubes or one 


of the G-M tubes was not functioning properly. 
These exceptions are indicated by the gaps in the 
graphs. An examination of these four curves 
clearly shows the marked irregularities in the 
north-west quadrant. The humps in the curves in 
this quadrant occur at the same angles in each 
curve. No such marked variations at consistent 
azimuth angles are shown in other quadrants. 
Curves A and C represent the data collected by 
opposite telescopes A and C during the same 
period, and curve E is the mean of A and C. 
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Similarly curve F is the mean of B and D and 
represents the data collected by two opposite 
telescopes B and D during the same period. In 
other words each point of curve F represents data 
taken about forty days after data taken for the 
corresponding point of curve E. In taking the 
means, it is assumed that the efficiency of the 
two opposite telescopes in each case is the same. 
By that the number of counts at each angle is 
doubled, and therefore the statistical errors are 
reduced. It should be noted that in the N-W 
quadrant all the points of curve F coincide with 
those of curve E within the probable errors ex- 
cept at 340°-350°. Even for these angles the 
general form of the curves is the same. The 
parallelism between curves E and F is clear 
enough, with the exception of regions between 
140°-170° and between 210°-250°. 

Curve G is the weighted mean of all the four 
telescopes. The means of all the ordinates of 
curve A, and similarly of curves B, C, and D, were 
taken separately. These means were weighted 
with reference to the mean of all the ordinates of 
the four curves (A, B, C, D) and it turned out to 
be the same as that of telescope A. That is, all 
the points of curves B, C, and D were multiplied 
by the ratio of their respective means to the 
general mean or that of curve A. It is again to be 
noted that in the N-W quadrant curve G has 
marked azimuthal irregularities. Two humps, one 
at 280° and the other at 310°, are prominent in 
this quadrant. The curve in the N-E quadrant 
varies considerably less. This is in agreement with 
Hutner’s conclusion that in the northern hemi- 
sphere, the curve between the N- and E-directions 
should be less irregular because of the supposed 
predominance of positive p articles. 


TABLE I, Counts per minute in N-W quadrant (N4) and in 
the S-E quadrant (Np). 


(~" -") 
a =2{ ——— 
Nat+Ns 





Angle NA Angle Nps 

0 0.272+0.005 180 0.298+0.006 —0.091+0.027 
350 0.263+0.005 170 0.280+0.005 —0.062+0.026 
340 0.285+0.005 160 0.273+0.005 +0.043+0.028 
330 0.262+0.005 150 0.267+0.005 —0.019+0.027 
320 0.274+0.005 140 0.267+0.005 +0.026+0.026 
310 0.280+0.004 130 0.290+0.004 —0.035+0.016 
300 0.301+0.004 120 0.278+0.004 +0.080+0.020 


290 0.301+0.004 110 0.272+0.004 
280 0.288+0.004 100 0.252+0.004 
270 + 0.276+0.004 90 0.256+0.004 


+0.101+0.020 
+0.133+0.016 
+0.075+0.016 
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Hutner’s results refer strictly to primaries be- 
fore they enter the atmosphere. In the present ex- 
periment the secondaries were filtered out by 10.2 
cm of lead, and therefore the telescopes recorded 
the passage of hard secondaries presumably 
having the same direction as the primaries which 
were responsible for producing them. 


THE EAST-WEST AND NORTH-SOUTH 
ASYMMETRIES 


Special cases of the ‘‘Azimuthal effect” are the 
E-W and N-S asymmetries. Tables I, II, and III 
give the experimental results on these asym- 
metries at Lahore (geomagnetic latitude 22° 
north) for a fixed zenith angle of 60°. Table I is 
the same as Table I already published,* except 
for the probable errors, which were not correctly 
calculated. Table II gives the results of the new 
experiment. (x) represents the counts per hour 
for each angle from N-E-S, while (y) gives the 
counts per hour from S-W-N. 


The asymmetry =2(x—y)/(x+y), 


is calculated at each 10° interval. 

A notable feature of these asymmetries is that 
the maximum E-W difference does not coincide 
with the magnetic E-W plane, but is found to 
occur in the 280°-100° plane. Each of the four 
telescopes gives the maximum E-W difference at 
the same angle (280°-100°). This becomes more 
significant when we look at Table I, which also 


TABLE II. Counts per hour from N-E-S(x) 


and from S-W-N(y). 
z—y 
(2) 
+y 











Angle x Angle y 
0 5.14+0.10 180 5.53+40.12 —0.073+0.029 
10 546+0.11 190 5.69+0.12 —0.041+0.029 
20 = 5.60+0.13 200 §=5.5240.12 +0.014+0.032 
30.» =5.44+0.13 210 =5.26+0.11 +0.021+0.032 
40 5.19+0.15 220. =5.12+0.13 +0.013+0.038 
50 §=4.96+0.13 230» =5.04+0.13 —0.016+0.037 
60 4.47+0.12 240 = 5.18+0.12 —0.147+0.035 
70 =4,98+0.13 250 5.1640.16 —0.035+0.040 
80 486+0.13 260 5.30+0.20 —0.087+0.046 
90 4.44+0.09 270» =5.40+0.12 —0.195+0.030 
100 4.62+0.09 280 5.80+0.09 —0.226+0.025 
110 . 4.83+0.11 290 5.61+0.11 —0.149+0.030 
120 4.92+0.12 300 4.99+0.10 —0.014+0.031 
130 4.7340.12 310 5.602+0.14  —0.167+0.035 
140 4.81+0.18 320 = 5.32+0.14 —0.100+0.045 
150 4.89+0.11 330 »=64.70+0.11 +0.040+0.032 
160 5,05+0.13 340 4.8840.10 +0.034+0.033 
5.10+0.12 350 5.29+0.16 — 0.037 +0.038 


170 

















350 y.jf. Ghee 
TABLE III. 
Angle A B Cc D E F G 
0 5.18=0.18 4.95=0.16 5.28 =0.22 5.15=0.22 5.23 =0.14 5.05 =0.14 5.14=0.10 
10 5.60 =0.25 4.97 =0.18 6.12 =0.30 5.17 =0.22 5.86 =0.20 5.07 =0.14 5.46=0.11 
20 5.55=0.23 5.22=0.18 5.98 =0.24 5.63 =0.33 5.76=0.17 5.42 =0.19 5.60 =0.13 
30 5.48 =0.23 4.84=0.21 6.17 =0.30 5.30 =0.32 5.82 =0.19 5.07 =0.19 5.44=0.13 
40 5.37 =0.22 4.53 =0.29 5.83 =0.25 5.60 =0.17 4.53 =0.29 5.19=0.15 
50 4.38 =0.20 4.46=0.29 5.57 =0.23 5.42 =0.32 4.98 =0.15 4.94=0.22 4.96=0.13 
60 4.11=0.18 4.25=0.28 5.18 =0.22 4.35 =0.29 4.64 =0.14 4.30 =0.20 4.47=0.12 
70 4.58 =0.30 4.74=0.21 5.29 =0.22 5.30 =0.32 4.94=0.18 5.02 =0.19 4.98 =0.13 
80 4.94=0.22 4.48 =0.17 5.13 =0.31 4.88 =0.30 5.04=0.19 4.68 =0.18 4.86=0.13 
90 4.17=0.15 4.28=0.20 5.15=0.18 4.17=0.20 4.66 =0.12 4.22=0.14 4.44=0.09 
100 4.40=0.14 4.59=0.16 5.00 = 0.16 4.50=0.22 4.70=0.11 4.54=0.14 4.62 =0.09 
110 4.68 =0.17 5.02 =0.22 4.88=0.21 4.70=0.21 4.78=0.14 4.86=0.15 4.83 =0.11 
120 5.28 =0.32 4.71=0.21 4.86=0.18 4.83=0.21 5.07 =0.18 4.77 =0.15 4.92=0.12 
130 4.87 =0.21 4.77 =0.22 4.54=0.21 4.77 =0.22 4.70=0.15 4.73=0.12 
140 4.75 =0.30 4.97 =0.31 4.42 =0.29 4.75=0.30 4.70=0.21 4.81=0.18 
150 5.17 =0.22 4.62 =0.21 5.32 =0.22 4.47 =0.18 5.24=0.22 4.54=0.14 4.89=0.11 
160 5.59=0.19 4.97 =0.31 5.52 =0.23 4.13=0.29 5.56=0.15 4.55=0.21 5.05 =0.13 
170 5.45=0.19 4.50=0.29 5.84=0.24 4.67 =0.21 5.64=0.15 4.58 =0.18 5.10=0.12 
180 5.52=0.18 5.16=0.19 5.71 =0.23 5.71=0.33 5.62 =0.15 5.44=0.19 5.53=0.12 
190 §.27 =0.14 5.88 =0.22 5.90 =0.24 5.88 = 0.22 5.58=0.14 5.69 =0.12 
200 5.52 =0.23 5.17 =0.18 5.93 =0.25 5.46=0.32 5.72 =0.17 5.32 =0.18 5.52=0.12 
210 5.19=0.22 5.32=0.18 5.46=0.23 5.07 =0.22 5.32 =0.16 5.20=0.14 5.26=0.11 
220 5.58 =0.23 4.44=0.27 5.50=0.23 4.97 =0.31 5.54=0.16 4.70=0.21 5.12=0.13 
230 5.04 =0.22 4.55=0.29 5.75 =0.23 4.84 =0.30 5.40=0.16 4.70=0.21 5.04=0.13 
240 5.10=0.22 5.45 =0.23 5.02 =0.19 5.10=0.31 5.06=0.15 5.28 =0.19 5.18=0.12 
250 5.42 =0.23 4.96 =0.22 4.96=0.22 5.42 =0.23 5.16=0.16 
260 5.42 =0.32 5.27 =0.24 5.27 =0.24 5.42 =0.32 5.30 =0.20 
270 5.26=0.19 5.26=0.23 5.40=0.17 5.66 = 0.33 5.33 =0.13 5.46=0.20 5.40=0.12 
280 5.73=0.15 5.79=0.21 5.85=0.19 5.80=0.18 5.79=0.12 5.80=0.14 5.80 =0.09 
290 5.65 =0.23 5.30=0.21 5.57 =0.19 5.88 =0.24 5.61=0.15 5.59=0.16 5.61=0.11 
300 5.13=0.22 5.13 =0.22 4.88 =0.18 4.79 =0.21 5.00 =0.14 4.96=0.15 4.99 =0.10 
310 5.23 =0.32 5.65 =0.23 6.00 = 0.34 5.54=0.23 5.62 =0.23 5.60 =0.16 5.60 =0.14 
320 5.21=0.32 5.68 =0.23 5.21=0.32 5.13=0.22 5.21=0.22 5.40=0.16 5.32 =0.14 
330 4.98 =0.22 4.83 =0.19 4.63 =0.30 4.32=0.20 4.80=0.18 4.58=0.14 4.70=0.11 
340 5.15=0.22 4.30=0.20 5.10=0.18 4.98 =0.22 5.12=0.14 4.64=0.15 4.88 =0.10 
350 5.58 =0.23 4.93 =0.30 5.46=0.32 5.52 =0.20 4.93 =0.30 5.29=0.16 


| 
| 
| 





shows that the maximum E-W difference occurs 
in 280°-100° plane. This point is being further 
investigated in this laboratory and it is hoped 
that the results will be reported in the near 
future. 

The present survey will be continued with a 
much larger apparatus, which is being con- 


structed for rapidly collecting new data. The 
author wishes to express his thanks to the Uni- 
versity of the Punjab for granting a fellowship 
for continuing these experiments at Lahore, to 
the University of Chicago for supplying the 
apparatus, and to Forman Christian College for 
the facilities of the Physics Department. 
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Demagnetizing Factors of the General Ellipsoid 


J. A. OsBorn 
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Charts and tables of the demagnetizing factors of prolate and oblate spheroids are readily 
available; however, demagnetizing factors of ellipsoids of three different axes are incompletely 
tabulated and laborious to calculate. This article presents charts and tables which make possible 
easy determination of the demagnetizing factor for any principal axis of an ellipsoid of any 
shape. Formulas for the demagnetizing factors of the general ellipsoid are included together 
with supplementary formulas which cover a large number of special cases. 


I, INTRODUCTION 


F a ferromagnetic body of irregular shape is 

brought into a uniform applied field Ho, the 
magnetizing force H inside the material differs in 
magnitude from the applied field and varies in 
direction throughout the body in an unknown 
manner. This is usually a great disadvantage in 
research on magnetic materials, because such 
research usually requires studies of the relation 
between the magnetization J and the magnetizing 
force. However, for homogeneous bodies whose 
surface is of the second degree,| H and J 
(after suitable magnetic treatment) are uniform 
throughout, though they are not necessarily the 
same in direction as Ho. The ellipsoid has the only 
surface of the second degree that is finite, so 
materials in this form are frequently used in 
precise investigations of magnetic materials. 

Inside any ellipsoid the component of H along 
any principal axis 7 is determined by the relation 


H;= (H)i—NiJFi, 1=X, Y, 2, (1.1) 


where NV; is a constant called the demagnetizing 
factor; it is determined by the ratios of the axes. 
Hence, if the magnetization components and the 
lengths of the axes of the ellipsoid are known, [7 
inside the body may be found.? 

The determination of N is rather involved, 
except for the special cases of ellipsoids of revolu- 
tion. Therefore, it is believed that calculation of 
N for various axial ratios of the general ellipsoid 


1]. C. Maxwell, Electricity and Magnetism (The Claren- 
don Press, Oxford, 1904), third edition, Vol. 2, pp. 66-70. 

* It is possible to find H at the center of a cylindrical bar 
if a kind of demagnetizing factor different from that for the 
ellipsoid is defined. The use of these so called ballistic 
demagnetizing factors is discussed by R. M. Bozorth and 
D. M. Chapin, J. App. Phys. 12, 320 (1942). 
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will serve a useful purpose. Naturally, it is 
impossible to include values of the demagnetizing 
factor for all axial ratios, so graphs of N have 
been made that will allow interpolation for those 
ratios not explicitly included. The following 
section gives the formulas from which the de- 
magnetizing factors have been calculated, to- 
gether with comments on the method by which 
they were obtained. 


II. FORMULAS FOR THE DEMAGNETIZING 
FACTORS OF THE GENERAL ELLIPSOID 


The equations given below are for the demag- 
netizing factors along the three axes of the general 
ellipsoid under the assumption that 


a>b>c>0, 


where a, 6, and ¢ are the ellipsoid semi-axes. The 
demagnetizing factors (corresponding to the semi- 
axes a, b, and c) will be labeled L, M, and N. 
Then, as indicated in Eq. (1.1), to find the 
component of the magnetizing force along any 
principal axis, the component of the applied field 
(Ho); along that axis must be considered together 
with the appropriate magnetization component 
J; and demagnetizing factor. The formulas for L, 
M, and N are 


cos ¢ cos 3 
L/4nr= 








[F(k, 3) —E(k, 8) ], (2.1) 


sin’ 3 sin? a 


cos ¢ cos 3 





- | ve, 3) 


sin’ 3 sin? a cos? a 


M/4nr= 





—cos? aF(k, 3) — 


sin? a sin 3 cos J 
| (2.2) 


COs ¢ 
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cos ¢ cos 3 fsin 3 cos ¢ 
N/4r= — E(k, | (2.3) 
sin* 3 cos? B cos 3 
where 
cos ’=c/a, (O<d8<-2/2), (2.4) 
cos g=b/a, (O<¢<7/2), (2.5) 
1—(b/a)*}! sing 
sin «-|— | = =k, 
1 — (c/a)? sin 3 
(O<a<n/2), (2.6) 


and F(k, #) and E(k, #) are elliptic integrals of 
the first and second kinds * k is the modulus, and 
#3 is the amplitude of these integrals. Equations 
(2.1)—(2.3) are given in forms that are convenient 
for calculation. Methods of obtaining general 
expressions for L, M, and N of the general 
ellipsoid are available in several places.‘ 
Maxwell has shown how the method of Poisson 
for getting the potential of a uniformly mag- 
netized body yields the demagnetizing factors for 
the general ellipsoid. The equations for the 
factors are not given in directly usable form how- 
ever. Kellogg‘ briefly indicates the method of 
obtaining L, M, and N and gives expressions 
similar to Eqs. (2.1)—(2.3). The most complete 
discussion of the solution for the demagnetizing 
factors of the general ellipsoid by use of the 
Poisson method is given by Poritsky.‘ Stratton‘ 
obtains expressions for L, 1, and N in the same 
form as found in Maxwell. 

A useful relationship that exists between the 
demagnetizing factors of an ellipsoid is 


L+M+N=4r. 


This may be proved by using the following ex- 
pression for the gravitational potential of an 
ellipsoid of uniform density p: 


V =(1/2)p(c—Lx?— My?— Nz’), 


(2.7) 


(2.8) 


where the coefficients L, M, and N represent the 


3See B. O. Pierce, A Short Table of Integrals, third 
edition (Ginn and Company, New York, 1929), p. 66 


et seq. 

‘6. Kellogg, Foundations of Potential Theory (Verlags- 
buchhandlung Julius Springer, Berlin, 1929), pp. 192-194, 
also p. 197, Exercise 7. H. Poritsky, Magnetization of an 
Ellipsoid (General Electric Company, Schenectady, New 
York). J. Stratton, Electromagnetic Theory (McGraw-Hill 
Book Company, Inc., New York, 1941), pp. 211-213. 

5 Poisson showed that if V is the gravitational potential 
of a body of uniform density p, the magnetic potential {, 
of the same body is —dV/dx, if it is uniformly magnetized 
in the x-direction so that J =p. 
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usual demagnetizing factors as in the Poisson 
method. Poisson’s operation 


eV #V a#V 
VV =—+—+—-=-—4rp, 


' (2.9) 
Ox® ay? a2" 











gives Eq. (2.7). 

If the ellipsoid is an ellipsoid of rotation, i.e., a 
prolate or oblate spheroid, Eqs. (2.1)—(2.3) reduce 
to more simple forms; the formulas for these 
and other limiting cases are given below: 


(a) If b=c 


| m 


m? —1L2(m?—1)! 


(prolate spheroid), 


L/4r= 








m+(m?—1)! 
x In (— — )-1} (2.10) 
m—(m*?—1)! 
where m=a/c. 
m 1 
M/42x = N/4x =———_ |m- _ 
2(m? —1) 2(m*—1)! 
m+ (m?—1)3 
xin ( )I (2.11) 
m —(m*—1)! 
If m>1 (very slender prolate spheroid), 
L/42r—=~(1/m?*)(In 2m—1), (2.12) 
M/4x=N/42r~3[1—(In 2m—1)/m?]. (2.13) 
(0) Ifa>b>c (very slender ellipsoid), 
L/4x = (bc/a*)(In 4a/(b+c)—1), (2.14) 
M /4x=c/(b+c) —4$(bc/a*) In 4a/(b+c) 
+bce(3b+c)/4a*(b+c), (2.15) 
N/42=b/(b+c)—}4(bc/a*) In [4a/(b+c) } 
+bce(b+3c)/4a2(b+c). (2.16) 
(c) Ifa=x,b>c (elliptic cylinder), 
M/4x=c/(b+c), (2.17) 
N/4x=6/(b+c). (2.18) 
(d) Ifa=b (oblate spheroid), 
1 
L/4e = M/42=————{ m?(m*?—1)=} 
2(m? —1) 
Xarc sin [(m?—1)*/m]—1}, (2.19) 








n 
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TABLE I. Demagnetizing factors of the general ellipsoid calculated from formulas of section I]. 
c/a b/a L/4r M/4r N/44 c/a b/a L/4r M/4r N/4e 
017452 99620 013471 .50000 1.00000 .23640 = .23640 52720 
98164 013719 98863 23555 —-.23885 52560 
95632 014165 95513 .23256 —-.24672 52072 
91357 014989 90139 22760 = .26012 51228 
40709 037634 83073 22043 = .27966 49991 
.29285 053078 .74825 21105 30560 48335 
19158 .080810 66144 19970 —-.33745 46285 
.071895 19374 58115 18772 —_.37190 44037 
055162 .23901 52213 17762 ~—- 40115 42123 
039018 .30809 50000 17356 41322 41322 
034900 89892 94467 64279 1.00000 = .27187 Ss .27187 45627 
.79891 94104 99111 .27097 .27403 45499 
69511 93617 96507 .26824 —-.28058 45118 
60246 93029 92374 26370 ~—-.29152 44478 
50091 92126 87037 25760 30650 43590 
40798 071847 90905 80971 24999. 32522 42479 
.29427 .099697 .88398 74825 .24148 34617 41235 
19386 14720 83986 69414 23330 = .36647 40023 
093838 .26746 .72453 65641 22717 ~—_ .38169 39114 
077961 .30596 64279 22486  —.38757 38757 
062880 35411 
1.00000 29 2 0660 
.087156 1.00000 .06154 .06154 .87692 — 99375 =. pote —— 
98492 06108 .06281 87611 97553 29388 —-.30303 40309 
.94013 .06051 .06620 .87329 94696 29048 31062 39890 
86727 05908 = .07222 -86870 91065 28595 32086 39319 
.76809 05683 08307 86010 87037 28064 33285 38651 
-64625 05347 -10028 84626 83073 .27513 .34529 37958 
.50565 04846 12912 -82242 .79697 .27018 .35644 37338 
35171 .04087 .18185 -77728 77413 26668 36464 36868 
.19376 .02898 .29661 .67441 76604 26542 36729 36729 
08716 01641 49180 49180 
17365 1.00000 11138) = .11138 = 77724 | "80003 = 1.00000 31380-31380 37240 
: . : : 99622 31352 31462 37186 
94157 10936 = .11903 77161 96824 "30083 32224 36780 
87037 10665 12955 "76380 (04696 30707 = 32803 36490 
77413 10245 14651 75104 92374 30396 33473 36131 
65641 09627 —-.17296 .73077 90139 "30093 «34121 35786 
52213 08734  .21496 69770 yr ’ 
~ + 88274 29825 - .34692 35484 
37895 07455 —-.28191 64354 87037 29651 35050 35200 
24371 05756 38888 55356 “86603 30584 «35208 35208 
17365 04582 47709 47709 a = . : 
25882 1.00000 15207 15207 69587 .93969 1.00000 32497 32497 35006 
98583 15136 = .15440 69424 -99823 32474 = 32548 34978 
94386 14927-16154 68919 ‘99314 32409-32673 34918 
87562 14553 17452 67995 98527 32306 32889 34805 
78392 13985 19464 66551 97553 32177 33153 34669 
67269 13167 22498 64335 .96507 32038 33435 34527 
54793 12029 26973 60998 95513 31902 33707 34391 
41970 10524 —-.33356 .56120 94695 31790 33947 34260 
30844 08817 41326 49857 ‘94157 3171634105 34183 
25882 07892 46054 46054 .93969 .31690 34155 34155 
34202 1.00000 18555 18555 62889 98481 1.00000 33119 33119 33762 
98660 18477 .18805 62718 99955 33120 33143 33737 
94695 18221 19597 62181 99824 33100 33180 33720 
88274 17755 .21045 61200 99622 33078 33230 33692 
.79697 17108 = .23128 59764 .99375 33044 33295 33661 
69414 16176 — .26200 57624 99111 33010 33364 33628 
58115 14940 —-.30440 54620 .98863 32977 = .33431 33590 
46933 13428 35885 .50687 98660 32950 — .33487 33563 
.37895 11924 —-.41586 46490 98527 32931 .33512 33558 
.34202 11218 = 44391 44391 98481 32925 33537 33537 
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TABLE II. Demagnetizing factors of the general ellipsoid for b/a equal to integral multiples of 0.1. 
b/a c/a L/4x M/4n N/4r b/a c/a L/4e M /4x N/4r 
1000 .0872 .0181 4624 .5194 .7000 .6428 .2341 3641 4018 
.0349 2545 .7359 .5000 .2050 .3226 4724 
0175 .1530 .3420 .1625 .2600 .5776 
.2588 .1338 .2168 .6494 
? 5 
2000 1736 0505 4400 5082 ‘Onn: 0851 0925 8598 
.0872 .0294 .2890 .6828 ‘0349 = nade ‘0364 
.0349 .1431 8439 re . 
0175 0785 8000 .7660 2706 3555 3739 
.6428 .2487 3284 4229 
3000 .2588 0866 4204 4927 -5000 2172 2887 4941 
1736 0653 3369 5990 3420 1713 2303 5984 
0872 0375 .2095 .7580 -2588 1408 -1908 6684 
0349 0977 8854 .1736 .1037 .1420 .7546 
0175 0518 .0872° .0576 .0798 .8630 
.0349 .9412 
.4000 .3420 .1230 4014 4760 .9000 .8660 .3008 3415 3577 
.2588 .1031 .3459 5522 .7660 .2845 .3240 3915 
.1736 .0768 .2692 6532 .6428 2611 .2981 4408 
.0872 0435 1613 .7940 .5000 2275 2605 5121 
.0349 .0727 .9083 .3420 .1788 .2063 .6148 
.0175 .0377 .2588 .1470 .1697 .6832 
.1736 .1077 .1250 .7670 
: .0872 .0597 .0691 8710 
.5000 .5000 1735 4132 4132 = 
3420 1387 3428 5186 0349 9445 
2588 1149 2920 5930 | 1.0000 0000 3333 3333 3333 
.1736 .0855 .2235 .6909 
0872 0481 1308 8217 .9848 3314 3312 3374 
a c 9397 .3249 .3250 3500 
.0349 .9202 
.8660 .3138 .3138 3724 
.7660 .2965 .2967 4066 
-6000 .5000 .1906 .3638 4457 .6428 .2718 .2718 4561 
.3420 .1516 .2962 .5521 .5000 .2364 .2364 5270 
.2588 .1253 .2501 .6253 .3420 .1856 .1858 .6288 
.1736 .0926 .1895 7173 .2588 .1520 .1520 .6960 
.0872 .0520 1080 8398 .1736 .1114 1117 7775 
.0349 9301 .0872 .0615 .0615 .8769 
N/4x = m?/(m?—1){1—1/(m?—1)! whose argument is 
Xarc sin [(m?—1)!/m]}. (2.20) e=(1—6?/a?)!. (2.26) 
If m>1 (very flat oblate spheroid), III. CONSTRUCTION AND USE OF TABLES 
AND GRAPHS 
L/4e=M/4ae—~(x/4m)(1—4/xm), (2.21 . 
/4n=M/4a=(x/4m)(1—4/4m), (2.21) Construction 
N/4n=1 — 4 /2m+2/m’. (2.22) The values of the demagnetizing factors for the 
(ce) Ifa>b>c (very flat ellipsoid), general ellipsoid given in lable I are determined 
by the quantities 3 and k (=sin a) used in Eqs. 
L/4 c — A-E 2.23) (2.1)—(2.3). These were chosen so that 8 and a@ 
/ — —’) os \é. covered the region from 0° to 70° in 10° steps, and 
from 70° to 90°, 3 varied in 5° and a@ in 10° steps. 
cE—(1—é)K Some gaps in the values of the demagnetizing 
M/4xr=- (2.24) : : - ; 
a e(1—e)! factors were revealed when the final graphs were 
plotted, as explained below; hence, additional 
N/te=1 cE 2.25 values were subsequently computed to make the 
/4a= all ~¢)! (2.25) graphs complete. The axial ratios of the ellipsoids 
were then determined by Eqs. (2.4)—(2.6). L/4x 
where K and E are complete elliptic integrals’ and N/4mr were determined directly from Eggs. 
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(2.1) and (2.3), and M/4x was then determined 
in most cases by means of Eq. (2.7). 

L/4x and N/4 in Table | are accurate to at 
least the fourth decimal place. M/4 is somewhat 
less accurate than the other factors because it 
involves a difference. 

The values of c/a and b/a obtained from 8 and a 
are such that interpolation for the demagnetizing 
factors of other ellipsoids is not easy. To reduce 
this difficulty, each factor was separately plotted 
on a large scale against b/a, with c/a as a pa- 
rameter. This gave a family of curves for each 
factor, L, M, and N. Each family of curves was 
then considered by itself. A value of b/a equal to 
an integral multiple of 0.1 was selected, and a 
value of the demagnetizing factor in question 
read from each curve. The constant value of b/a 
was now taken as a parameter and a new curve 
plotted with the factors (read from the above 
family of curves) as ordinates and the values of 


c/a (the parameter of the above family of curves) 
as abscissas. This process was repeated for 
different b/a values until a new family of curves 
was obtained with the parameter }/a equal to 
integral multiples of 0.1 from 0.1 to 1. These 
curves (shown in Figs. 1-3) warrant graphical 
interpolation to three decimal places. Table |! 
contains the values of L/4x, M/4x, and N/4x for 
various c/a (with b/a as a parameter) from which 
Figs. 1-3 were plotted. These values are accurate 
to three decimal places and are probably in error 
several units in the fourth place. 


Use of Graphs 


To use the graphs, any value of c/a may-be 
chosen within the limitations 


a>b>c>0 


and b/a may be found by interpolation. It is 
most important to remember that: L corresponds 
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to a, the longest semi-axis; M corresponds to b, the along length 3 is wanted. Here 
intermediate semi-axis; and N corresponds to c, the e=3: b=2?: and c@1. 
shortest semi-axis. Two examples are given below 
illustrating the procedure for use of the graphs. 
Examples: (1) Suppose the semi-axes of an 
ellipsoid are 4, 3, and 2 and it is desired to find 
the demagnetizing factor along the axis of 
length 3. Then If the magnetization curve is linear, and if the 
material either is isotropic or has its principal 
axes of magnetic anisotropy along the principal 
The factor is then read from the curves for M/4r axes of the ellipsoid, then 

with axial ratios c/a=0.500 and b/a=0.750. Leaiieln ~tilidie. (4.1) 
M /4x =0.306 when read from the curve. (2) Sup- 
pose the axis lengths are 3, 2, and 1 and the factor 


In this case the factor is read from the curves for 
L/4x with axial ratios c/a =0.333 and 6/a =0.667. 
L/4r =0.156 when read from the curve. 


IV. APPLICATION TO SPECIAL MATERIALS 


az4:. b=3: and c=z2. 


where the susceptibility x; and the permeability 
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ui, Corresponding to axis i, are constant. (For By use of these formulas, the magnitude and 
isotropic materials the three x;’s are equal.) direction of the magnetizing force and mag- 
Therefore netization may be found if the magnitude and 


H;=(He):—NiJi=(He)i—Naxil:, (4.2) direction of the applied field are given. 
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A simplified discussion is given of the problem of introducing linearly accelerated axes by a 
method applicable alike to classical and relativistic mechanics. The mathematical method is 
based on the Lie theory of the 4-dimensional conformal group which is presented in such a way 
that the transformations leading to accelerated axes are exhibited explicitly. The classical 
theory is shown to be a degenerate case of the relativistic theory, so that the transformations to 
accelerated axes have a geometrical interpretation in terms of inversive transformations in both 


theories. 





1, INTRODUCTION 


ROM time to time, modifications of the special 
theory of relativity to include transforma- 
tions between relatively accelerated Euclidean 
coordinate systems have been discussed in the 
literature. The first steps in this direction seem 
to have been taken by Einstein' and Born.? More 
recently, the question has been investigated again 
by Milne® in connection with his theory of the 
expanding universe, while Paget has proposed a 
detailed theory of the kinematical side of the 
problem. The work of Milne and Page has been 
subjected to critical discussion by Bourgin,® and 
more particularly by Robertson® who has given a 
very penetrating treatment of the general mathe- 
matical theory, with particular reference to its 
implications in cosmological theories. Finally, 
Engstrom and Zorn’ have pointed out inde- 
pendently the connection between Page’s theory 
and the 4-dimensional conformal transformation 
group. Unfortunately, the widely divergent phys- 
ical and mathematical viewpoints of the various 
authors make it difficult to arrive at a clear view 
of the physical implications of their discussions. 
The present work originated in an independent 
attempt to develop a relativistic equivalent of the 


! A. Einstein, Jahrbuch der Radioaktivitat 4, 411 (1907). 

2 M. Born, Ann. d. Physik 30, 1 (1909). A brief discussion 
is also given by H. Bateman, Proc. Lond. Math. Soc. [2] 8, 
223 (1910). 

3E. A. Milne, Relativity, Gravitation and World-Structure 
(Oxford University Press, New York, 1935). 

*L. Page, Phys. Rev. 49, 254 (1936). 

° D. Bourgin, Phys. Rev. 50, 864 (1936). 

*H. P. Robertson, Phys. Rev. 49, 755 (1936); Zeits. f. 
Astrophys. 7, 153 (1933); Astrophys. J. 82, 284 (1935), 83, 
187 (1936), 83, 257 (1936). 

7H. T. Engstrom and M. Zorn, Phys. Rev. 49, 701 
(1936). 


simple process in Newtonian mechanics whereby 
one introduces a system of coordinates moving 
with an accelerated particle. In this paper we 
shall give a direct solution of this problem by the 
simplest applicable means and in such a manner 
as to bring out the complete analogy between the 
classical and relativistic theories. 

While the present discussioh is non-quantum 
mechanical in character, it may be noted that the 
establishment of wave equations invariant under 
the conformal group has been discussed in the 
literature.® 

2. STATEMENT OF THE PROBLEM 


The starting point of our discussion will be the 
replacement of the ultimate problem of estab- 
lishing transformations between accelerated co- 
ordinate systems by another which at first sight 
may perhaps appear to have only a trivial 
significance. However, its solution gives the key 
to the more general question. - 

Apart from the trivial case of uniform motion 
under no forces, the simplest situation in me- 
chanical theory for which a clear-cut formulation 
can be given without involving difficulties in the 
discussion of force fields is that of the uniformly 
accelerated motion of a particle. This can be 
characterized in both classical and relativistic 
theories as a motion in which the acceleration is 
constant when measured in a system of coordi- 
nates in which the particle is instantaneously at 
rest. Such a set of coordinates is called a rest- 
system; a set of coordinates in which the particle 
is not only at rest but also has no acceleration 


will be designated as a proper-system. 


§P. A. M. Dirac, Ann. Math. [2] 37, 429 (1936); H. J. 
Bhabha, Proc. Camb. Phil. Soc. 32, 622 (1936). 
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Since a particle at rest is but a special case of 
one in uniformly accelerated motion, the trans- 
formation from a rest-system to a proper-system 
has the property that it transforms one type of 
uniformly accelerated motion into another. This 
suggests the following formulation of the problem: 
To determine those transformations for which the 
uniformly accelerated motion of a particle is trans- 
formed into another of the same type. 

In this form the problem can be attacked in the 
same manner in both classical and relativistic 
theories, and it has the further advantage of 
suggesting of itself the most direct method of 
solution. In the discussion of the coordinate 
transformations which arise, we shall proceed on 
the consideration that each such set of coordi- 
nates may be interpreted as associated with an 
“‘observer’’ to whom the motion of the particle as 
described by his coordinates is its true motion. 


3. THE CLASSICAL ONE-DIMENSIONAL PROBLEM 


In order to clarify the essentials of the method, 
we start with the simplest possible classical case. 
We consider a particle which in some appropriate 
initial coordinate system moves along a straight 
line with a constant acceleration g. Using x as its 
coordinate, the differential equation of its motion 
is 


d*x/dt? =g. (1) 


The differential equation characterizing all 
uniformly accelerated motions is therefore 


d*x/d*®=0. 
With the notation 
v=dx/dt, a=dv/dt, b=da/dt, 
this characteristic differential equation becomes 
F.(x, t, v, a, b) =b=0. (2) 


Our problem is now reduced to that of seeking 
all transformations of coordinates which leave 
this equation invariant. In each such new system, 
the particle will appear to move with constant 
acceleration, but the magnitude of the accelera- 
tion will vary with the system. | 

A study of the transformation properties of 
I:q. (2) shows that the only point-transformations 
which leave it invariant are those generated by 
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the following set of operators :* 


X,=—90:, X2=— 91, 

X;= —/0,, X,=—}9., (3) 

X5=—xtd,—}300,, Xg=xd,+10:, 
X7=xd,—(0:, 


where, for simplification, 0, =0/0x, 0:=0/dt. 
For convenience we give also the formulas for 

the general infinitesimal transformation, ex- 

tended three times to cover the order of Eq. (2): 


x’ =x—a,—la3— }fay—xlagst+Kxagt+Xaz, 

t! =t—a2—}Pagt+tasg—taz, 

v’ =v—a3—lag—xag+2vaz, (4) 
a’ =a—a4— (v—la)as—das+3aa;, 

b’ = b+2tba;— 2bag+4baz. 


Since these are extended point-transformations, 
the formulas for v, a, and 6 are consequences of 
those for x and ¢. 

The invariance of Eq. (2) follows at once from 
the fact that the quantity (b’—0) vanishes with 
b; i.e., vanishes when Eq. (2) is satisfied. 

The requirements of classical theory now lead 
us to anticipate the possibility of finding a 
4-parameter subgroup of transformations which 
will include one of the type sought to an ac- 
celerated coordinate system. The four parameters 
would be used, (a) to permit shifting the origin in 
the (x, ¢)-plane in order to give the particle the 
position x =0 at a given instant ¢=0 (2 parame- 
ters), (b) to transform to a rest-system with the 
particle at the origin (1 parameter), and finally, 
(c) to transform to the proper-system with the 
particle still at the origin (1 parameter). It is to 
be expected also that this can be achieved using 
only transformations for which time intervals are 
invariant. 

A study of Eqs. (4) shows that this can just be 
accomplished by using the four transformations 
generated by X1, X2, X3, X4. It follows from the 
results to be obtained in the next section that 
these transformations form a 4-parameter group. 
The finite transformations of this group can be 


* Throughout the paper the operators have been written 
in such a manner as to facilitate the physical interpretation 
of the resulting transformations. 
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calculated readily'® and reduced to the form 
x’ = (x—xX9) —v0(t—to) —Fao(t—to)?, t’ =t—to, 
with the notation 
Xo = ait (Saras+ayay”)/6; to=a, 
Vo=astaras/2; do=ay. 


In this form it is immediately recognizable as a 
change to moving, accelerated coordinates with 
superposed changes of origin for x and ¢t. By a 
proper choice of x» and ¢o, we bring the particle 
to the origin in the (x, ¢)-plane ; by equating vo to 
the velocity of the particle, we obtain a rest- 
system ; and by equating a» to the acceleration of 
the particle in the rest-system, we reduce it to a 
proper-system. For general values of do, we 
obtain a 1-parameter family of accelerated axis 
systems. 


4. THE RELATIVISTIC ONE-DIMENSIONAL 
PROBLEM 


In the special theory of relativity, the equation 
of uniformly accelerated motion, replacing 
Eq. (1), is 

ea/(?—v*)i=g, (5) 
where g is the constant acceleration as measured 
in a rest-system. By a further differentiation we 
find the differential equation characteristic of all 
uniformly accelerated motions to be 


F,(x, t, v, a, b) =b+3va?/(c?—v?)=0. (6) 


The point-transformations which leave this 
equation invariant form a 6-parameter group 
built on the operators :" 


X:=-0,, X2=—91, 

X3= —t-0,—(x/c*)-d:, 

X4= — (CP +2") /2C-0,—(xt/c?)-d., (7) 
X;,= —xt-0,—(CP +x") /2c?-d:, 
X¢=xd,+10:. 


10 The necessary theory of continuous groups can be 
found in the book of J. E. Campbell, Theory of Continuous 
Groups (Oxford University Press, New York, 1903) or that 
of G. Kowalewski, Einftihrung in die Theorie der Kontinuier- 
lichen Gruppen (Akademische Verlagsgesellschaft, Leipzig, 
1931). The method of finding the finite transformations of 
a group which is defined in terms of its infinitesimal ele- 
ments is discussed in Campbell, p. 47. 

The full set of symmetry transformations of both 
Eqs. (2) and (6) form a 10-parameter family, of which the 
remaining ones are contact transformations, but not point- 
transformations. 
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TABLE I. Commutator table of the symmetry operators 
of relativistic one-dimensional uniformly accelerated 


motion. 






































xX 0 

X 0 0 

Xs i Xai 0 

X, |-X/el xX; | Xs/e re 

zm. i Rs | Me x, 0 ae 

X. | -%¥:1-X,| © a : Xs 0 
x my xX, Re ks 


For convenience of reference we tabulate the 
general infinitesimal transformations to the third 
orders of derivatives: . 


x’ =x — a) — a3t —ay(CP +x?) /2c? — asgxtt+ agx, 


‘= t— a2 — asx /c?—ayxt/C 
ne a;(c?t?+x?) /2c?+ asl, 


v’ =v—a3(c?—v")/c?—agt(C —v’)/c? 
—asx(c?—v*) /c*, 
, | 9 9 9\ /.9 (8) 
a’ =a+a330a/c?+a,4(xa+ 3tva+v* —c*) /c* 


+as[ cta+3xva — cv +" ]/c? — aca, 
b’ =b+.a;(4vb+ 3a?) /c 
+ a,(2xb + 4tvb + 6va + 3ta*) /c? 
+ a;5(2c*th + 4xvb + 6v*a + 3xa*) /c? — ag2b. 


It is readily verified that this transformation 
leaves Eq. (6) invariant. 

Table I is the commutator table for the set of 
operators of Eq. (7). 

Each entry in Table | gives the commutator 
of the operator at the side of its row with that at 
the bottom of its column. The operators X,, X 2, 
X;3 generate the 3-parameter subgroup of transla- 
tion and velocity transformations, while X,4, X5, 
Xe generate a second 3-parameter subgroup. 
Neither of these subgroups is normal within the 
whole group. 

In order to obtain the operators of Section 3, 
we need only take the limit c— © in the operators 
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of Eq. (7), in the general infinitesimal trans- 
formations of Eq. (8), and in Table I.” 

To develop the ‘‘acceleration transformations” 
we now suppose that the starting coordinate 
system has been chosen as a rest-system in which 
the particle is brought to the origin in the 
(x, t)-plane. Equations (8) show that the sub- 
group of transformations generated by X4, X5, X¢ 
has the characteristic that the particle is left at 
the origin at rest. 

The 1-parameter subgroup of transformations 
generated by X, alone gives a family in which the 
particle is left at the origin at rest, but for which 
its acceleration becomes (g—a,). By equating a, 
to g, the proper-system of the particle is found. 
The finite transformations of this family are 
easily computed and found to give exactly the 
result found by Page‘ in his Eqs. (27) and (28), 
with a,=¢." 


5. GEOMETRICAL INTERPRETATION 


Each of the transformations used can be given 
a geometrical interpretation as an operation in the 
(x, t)-plane: (a) X; and X_ generate translations 
of the origin, (b) X3 has the usual interpretation 
given in the special theory of relativity as a real 
rotation in the complex (x, ict)-plane or as an 
imaginary rotation in the real plane, (c) in the 
real plane we define the transformation 





R°(x — Xo) 
T (xo, to): x’ —xo= — 
; (x—x9)?—C(t+ty)? 
Rit —to) 





t! —t5=—— — 


(x —x9)? —c?(t—to)? 


2 The operator X; of Section 3 appears in the relativistic 
theory as a degenerate case of a contact transformation 
which becomes a point-transformation on taking the limit 
c—>e, The group used here is isomorphic to the symmetry 
group of the circles in the complex (x, ict)-plane. The 

ssage from the relativistic to the classical theory appears 
in Lie’s geometrical theory as the transition from the 
symmetry group of circles in the complex plane to that of 
the family of parabolas in the real plane forming the 
classical trajectories of possible uniformly accelerated 
particles. This geometrical theory is discussed in the books 
of Campbell and of Kowalewski mentioned in reference 10. 

‘8 The finite transformations of the 2-parameter subgroup 
based on X, and X; have been computed. The equations 
show three singular lines, of which two have the equations 


(x+&ct)(aytasc)+2c?=0 


and are the generalizations of those noted by Page. The 
third is the line 


asx —agl =0 
which transforms into itself according to the relation 
(asx’ — agt’) « (x? — CF) = (asx —ard) - (x? —c*t"2). 
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as an inversion in the hyperbola 
(x—Xo)? —Ci(t —to)? = R?. 


If we let S(xo, to) be a translation of the origin 
through the distances x» and t») along the axes, 
then the infinitesimal transformations ayX, and 
asX,5 are given by the following sequences 


S(— R’a,/2c*, 0) 
-T(— R’a,/2c, 0)-T(0, 0), 
S(0, Rtas/2c?)-T(0, R°as/2c*) - T(0, 0), 


A4A4: 


A5A 5: 


the order of the transformations being read from 
right to left. In the complex plane, these become 
combinations of inversions in circles combined 
with translations,“ and (d) the transformations 
generated by Xz are uniform stretchings and 
contractions of the (x, ¢)-plane. 

From this we see that, even in the classical 
theory, the transformations to accelerated co- 
ordinate systems have a geometrical interpreta- 
tion in terms of inversions. This result carries 
over unchanged into the 3-dimensional theory. 


6. THREE-DIMENSIONAL MOTION 


The characteristic differential equation for 
3-dimensional uniformly accelerated motion is so 
complex that a direct study of its invariance 
would be very laborious. It is not difficult to set 
up the generalization of the operator system of 
Section 4 by direct evaluation, but it can be put 
on a more elegant basis by recognizing at once the 
connection with Lie’s theory of the group of 
conformal transformations in 4 dimensions. The 
physical basis of the connection is that since a 
particle traveling with the velocity of light repre- 
sents a special case of uniform acceleration in 
which further acceleration is impossible, the 
family of paths of light rays must therefore be 
invariant under the transformations. Since the 
differential equation of a light surface in the 
starting system is 


(dx)?+ (dy)?+ (dz)? —c(dt)? =0, 


this relation must also be satisfied after the 
transformation, which implies that 


(dr’)?—c2(dt’)? = p[ (dr)? —2(dt)*}, 


A discussion of this 2-dimensional geometrical theory 
can be found in Chap. 19 of Campbell, reference 10. 
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TABLE II. Commutator table of the conformal group. 


















































































































































Xi 0 
a See 
N: 0 | 0 
\ 0 | 0 0 
| i 
N« 0 | 0 0 0 
—| —— 7 
X 0 | X3 —X 0 0 
|__|} 
Xé —X; 0 xX 0 —X; | 0 
mk es Me o— 
| 
X X | —Xi 0 0 X¢ —X 0 
sain i ——| | 
X Xa/c? 0 0 X 0 X10 —Xs 0 
X» 0 X4/c? 0 Xe —Xw 0 Xs X: c 0 
| 
X10 0 0 X4/c? Xs; Xo —Xs 0 —Xe/c? X/¢ 0 
i ccdtal — = . | 
Xu —Xi;/c2| X7/c? | —X6/c? Xs 0 Xis —Xi Xu/c? 0 0 0 
ail C4 Stel a : Bata 
Xu —Xi/c? | —Xis/c?|. X5/c? Xs —Xi3 0 Xu 0 Xu/c? 0 0 0 
Xia Xe/c? | —X3/c? | —Xus/c*] Xw Xn —Xu 0 0 0 Xu/e 0 0 0 
Xu Xs Xo Xi —X, 0 0 0 Xu X) Xs 0 0 0 0 
X —X: | -X: | —Xs | -X. 0 0 0 o | oO 0 Xu Xu Xu | Xu 0 | 
xX 1 X 2 x 3 Xs Xs Xe X; Xs Xs Xvw Xu Xn Xis Xu Xi , 


where p is an unspecified function of x, y, 2, ¢. 
This relation may be used to characterize the 
group of transformations.” !5 

We proceed at once to give the set of 15 
operators which form the basis of the group of 
transformations: 


X,:=-0,, X2:=-0,, X3=-0, Xy=-O%, 
X5=—yd.+20,, X¢=—20,+x0,, 
X;=—xd,+y0., Xs= —ld,—(x/c*)d;, 
X9= —td,—(y/c?)d., Xyw=—td,—(2/c*)d,, 
Xn=— (CP +x? — y? —2*) /2c?-d, 

— (xy/c?)d, — (xz/c?)d, — (xt/c*)d,, 
X12= — (xy/c?)d.— (CP —x?+y? — 2") /2c?-d, 

— (yz/c*)d,—(yt/c?)d:, 


%S. Lie, Theorie der Transformationsgruppen (Teubner, 
Leipzig, 1930). A very interesting, but difficult, discussion 
is given in a classic paper by H. Bateman, Proc. Lond. 
Math. Soc. [2] 8, 223 (1910). 


X13= — (x2/c*)d,—(y2/c*)d, 
— (cf —x? —y?+2*) /2c?-d,—(2t/c*)d:, 


X 4= —xtd, — yld, — 20, — (CP +2°+y°+2*)/2c?-d,, 
X 15=xd,+ydy+20.+10:. 


The commutator relations of these operators 
are given in Table II. 

X1, X2, X3, X4 generate the translations. X5, 
Xs, Xz, generate the 3-dimensional rotations, 
while Xs, X9, Xi generate the Lorentz trans- 
formations. X11, X12, X13 form the 3-dimensional 
acceleration sub-set, of which X ,4 gives the time- 
component. X,,; again gives the homogeneous 
dilatations of the whole space. 

The procedure for studying the transformation 
to accelerated axes proceeds by the same princi- 
ples as in Section 4. By extending the trans- 
formations as far as the acceleration terms, one 
finds that if the particle is brought to the origin 
of coordinates in a rest-system, then it will be 
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left in this condition under all transformations of 
the acceleration subgroup; i.e., under the 3- 
parameter family of transformations generated 
by X11, X12, X13." 

On orienting the x-axis along the direction of 
the acceleration of the particle, the 1-parameter 
family of transformations which include its 
proper-system can be found. The transformation 
equations are: 


x’ =[x+an(r?—CF)/2c?)}- (t/t), 

y=y- (t/t), 2 =2-(t'/b), 

U=t/ {1+ (an/c?)(x+(an/4e)(r?—Cr) }}, 
r=x?+¥?+2". 


The passage to the classical theory is readily 
made, as in Section 4, merely by taking the limit 
c— in all of the equations for the operators as 
well as in the commutator table. The general 
infinitesimal transformation may then be written 
in 3-dimensional vector form as 


ér=r' —r= — bs— bw Xr — tov — 35a —rlayyt+rais, 
b= —t= —a4—Flanttars, 
with the abbreviations 

5S = (a,i+a2j+ask), 

dw = (asit+acsj+azk), 

dv = (agit+agj+aiok), 

6a = (a@),i+aj2j +a3k). 


The operators X ,4 and X,,; can be dropped, and 
the transformations to accelerated axes based on 
‘6 The particle at the origin in a rest system is left at that 


position at rest under the 5-parameter group generated by 
X ly X 12, X 13, X 14, X 15- 








the remaining 13-parameter group. Time inter- 
vals are then invariant. 


7. ELECTRODYNAMIC CONSIDERATIONS 


The detailed proof of the invariance of the 
electromagnetic field equations and the deriva- 
tion of the transformation relations of the field 
vectors and potentials have been given by 
Cunningham" and Bateman.'® However, this 
only assures us that electromagnetic phenomena 
as observed by any of the family of observers 
related by transformations of the conformal 
group will appear to obey the field equations. 
Since the transformations are in general non- 
linear, plane waves do not transform into plane 
waves, and the existence of singular regions in the 
finite transformations makes it unclear to what 
extent retarded and advanced potentials may be 
intermixed. 

We shall note here in conclusion that if one 
assumes that in its proper-system a point charge 
exhibits its usual Coulomb field, then Bateman’s 
formulas can be used to find its field in other sets 
of axes. A first-order calculation of the field due 
to a charge e at the origin of coordinates and 
moving with infinitesimal velocity 6v and ac- 
celeration 6a yields the result 


E=er/r’—e(r- da)r/2r°c* —eda/2rc’, 
H= —e(rX dv) /r*c. 


which is in agreement with the usual formula to 
this order of approximation. 

It is of interest to note also that if the po- 
tentials are computed, they contain contributions 
in ai, and as, even to first orders, but on passing 
to the field vectors these terms disappear, so that 
they introduce a type of gauge transformation. 


17 E, Cunningham, Proc. Lond. Math. Soc. [2] 8, 77 
(1910). 
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The Magnetic Field Inside a Ferromagnet 


G. H. WANNIER 
University of Iowa, Iowa City, lowa 
April 13, 1945 


T is generally agreed that for a fast charged particle (in 

the following, referred to as the “test charge’’) trav- 
ersing a magnetized iron bar, the average magnetic field b 
is equal to the induction B.“? However, doubts as to the 
correctness of this statement have been expressed on 
classical grounds* because the dominating contribution to 
this average is made by the enormous field which exists 
inside the spinning electron. This same state of affairs 
prevails in quantum theory. It expresses itself in the 
formula 


b=h+4xpM. (1) 


Here b is the average field, M is the magnetization, and h 
is the field which would exist if the electrons were true 
magnetic ‘‘di-poles.”” p is a numerical factor which is equal 
to the relative probability of coincidence of test charge and 
electron, as compared to randomness. 

The formula can be derived without specifying the struc- 
ture of the wave function, except for the following as- 
sumptions: 

(a) that the magnetic interaction is sufficiently small 
to be treated as a first-order. perturbation, 

(b) that the magnetic field is owing entirely to the 
electronic spin, 

(c) that the ferromagnetic electrons move in orbits inde- 
pendent of each other, ‘ 

(d) that it is sufficiently accurate to solve the Dirac 
equation in the Schroedinger type approximation of 
Darwin-Pauli, and 

(e) that the test charge is much heavier than an electron. 

Formula (1) becomes standard when the coincidence 
probability p equals 1, that is for plane waves. However, 
the term 47M _ is owing entirely to head-on collisions. This 
gives rise to two objections:* 

(1) Head-on collisions may not be frequent enough to 
average out for each particle. This objection is actually 
unfounded as such a collision occurs in a crystal about once 
every hundred atoms. 

(2) The coincidence probability p of the test charge and 
electron may be different from unity. 

The latter proposition is certainly true to some extent, 
partly because of the crystalline field. But the primary 
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TABLE I. 
p 
v/|Z| +charge —charge 
4-109 5.58 0.021 
$-10° 2.97 0.184 
1-10° 1.85 0.461 
2 -10° 1.39 0.692 
4-10° 1.18 0.837 








reason is the direct attraction or repulsion between electron 
and test charge upon close approach. Barring the case of 
uncharged particles such as neutrons, or the existence of 
short range forces yet unknown, this direct interaction is 
primarily the Coulomb interaction between electric charges. 
For this case, the numerical value of p is known.‘ It equals 


p=x/(1—e7*), (2) 
where 

x =49°Ze?/hy = +1.39 X 10°| Z| /2. (3) 
Here 7 is the velocity of the test charge in cm/sec. and Ze is 
the charge of the test charge. The sign of x depends on the 
sign of this charge. A few values of p are listed in Table I. 
Formula (1) suggests that if b is different from B, or h 
different from H, it will be primarily caused by short range 
Coulomb interaction between the test charge and the 
ferromagnetic electron. Such a calculation can be handled 
as a collision problem for two particles which are free 
otherwise. The necessary integrations can then be per- 

formed exactly and give 
b=H+2x(p+1)M, (4a) 


or 


b=B+2x(p—1)M, (4b) 


where p is again the parameter used in (1) and (2). This 
gives a field larger than B for positive charges and smaller 
than B for negative charges. 

From an experimental point of view, the most hopeful 
feature for a check is the assymmetry of the equations for 
positive and negative charges. Positive particles should be 
deflected more under otherwise similar circumstances, 
provided their speed is not much greater than 510° 
cm/sec. It should be added in conclusion that because of 
assumption (d), all formulas are subject to relativistic 
modifications which are being investigated by the author. 

1 C. F. v. Weizsaicker, Ann. d. Physik 17, 869 (1933). 

2? F. Rasetti, Phys. Rev. 66, 1 (1944). 

3W. F. G. Swann, Phys. Rev. 49, 574 (1936). 


4 Mott and Massey, Theory of Atomic Collisions (Oxford University 
Press, New York), p. 36, formula (17). 





Development of Electromagnetic Theory for 
- Non-Homogeneous Spaces—A Correction 


BENJAMIN LIEBOWITZ 
New York, New York 
May 5, 1945 


M* paper under the above title! purported to effect a 
separation of the components which define an 
electromagnetic field for the general case in which ¢, the 
specific inductive capacity of the space, is an arbitrary 
twice-differentiable function of the coordinates. The 
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method consists in introducing a Riemannian line element 
(Eq. (2) of the paper) with h?=1/e, and also in introducing 
an auxiliary scalar function @ for the purpose of obtaining 
the desired generality. This function Q is related to the h 
of the line element by Eq. (11). It has been called to my 
attention, however, that terms in the second derivatives of 
© were omitted in my Eq. (7); these terms cancel the Q-term 
in Eqs. (11) so that that function is not determined. When 
the term in @ is dropped in Eqs. (11), the theory remains 
valid only if the h of the liné element is restricted to the 
form of a product of functions, viz. 


h= A,(u;)H2(ue)H3(us), (A) 


where the H’s are arbitrary functions of their respective 
arguments with required continuity. Thus, equations such 
as (12), (18), (20), and (22) remain valid if 2 is put equal 
to unity and A is restricted in accordance with Eq. (A) 
above. However, wave equations in which h is put equal to 
a function of r (r? = u,?+u4*+-u;*) are valid only for the very 
special case in which the function is exp kr®. All my dis- 
cussion based on putting 4 equal to an arbitrary function 
of r is no longer justified. 

Although Eq. (A) still allows considerable latitude in the 
choice of h, nevertheless, the restriction it imposes is too 
severe from the standpoint of atomic theory, for which the 
work was originally intended. Attempts to restore greater 
generality by introducing a broader line element, viz., 


ds* = (fdu;? +fe2dus? + fsdu5*) /e, 


do not appear to help the situation materially, because, in 
order to effect the separation, it appears necessary to 
restrict the coefficients so that 
Si=fi(us), €=gi(ui)ge(us)gs(us), 

where each g and each f is an arbitrary function of the 
indicated coordinate only. I hope, however, to be able to 
report better success in the special case of a spherical 
coordinate system with ¢ a function of r only. 

I am greatly indebted to Dr. T. Pearcey of England for 
calling this error to my attention. 


Phys. Rev. 64, 294 (1943). 





On the Evidence for the Existence of 
Neutral Mesotrons 


S. C. SIRKAR AND P. K. BHATTACHARYYA 
Department of Physics, University College of Science, Calcutta, India 
June 2, 1945 


T has been observed by Groetzinger and others! that 

when 10.5 Mev deuterons strike targets of Mn, W, P, 
Sb, or NaCl two-fold coincidences occur in a counter tele- 
scope placed outside the cyclotron at a distance of about 
four meters from the target in the same line as that of the 
deuteron beam, even when the cyclotron is shielded by 
water tanks 120 cm thick and the counter telescope is 
shielded by lead sheets 19 cm thick. The walls of the 
counters were such that only 2 Mev electrons could pass 
through both the counters. From results of investigations 
on the influence of the thickness of the lead shield and 
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water shield on the number of coincidences observed, they 
have concluded that the electrons responsible for the ob- 
sérved coincidences are produced by the decay of some 
neutral particles coming from the targets bombarded by 
10.5 Mev deuterons, and these particles are probably low 
mass neutral mesotrons. 

Evidence for the presence of such unstable neutral 
particles in cosmic rays was found in 1942 in some in- 
vestigations carried out in this department? with an 
anti-coincidence counter telescope placed first with its axis 
vertical, and next inclined to the vertical at an angle of 30°. 
It was observed that some neutral particles, passing through 
the anticounter without discharging it, produce charged 
mesotrons in a lead absorber 5 cm thick and these mesotrons 
produce a three-fold coincidence in a counter telescope. The 
intensity of such neutral particles was found to be about 
1.4 percent of the total mesotron intensity in the laboratory, 
the effect due to side showers being eliminated completely 
The intensity was also found to diminish in the inclined 
direction almost in the same ratio as that of the charged 
mesotron. The apparatus used is shown diagrammatically 
in Fig. 1. The absorption coefficient of these particles in 











Fic. 1. Arrangement of apparatus in the counter telescope. 


lead was found to be less than that for charged mesotrons. 
From these results it was concluded that these neutral 
particles are unstable, and they have a life of the same 
order as that of the charged mesotron. It can be easily seen 
now that the particles observed by Groetzinger and others 
to be produced artificially by bombardment of high energy 
deuterons are identical with those found to be present in 
cosmic rays. The investigations of Groetzinger and others' 
confirm the inference made from results of cosmic-ra\ 
investigations mentioned above that these particles are 
unstable and give rise to high energy electrons after 
spontaneous disintegration. 

The writers of the note are indebted to Professor M. N. 
Saha, F.R.S. for kindly drawing their attention to the note 
by Groetzinger and others. 

1G. Groetzinger, P. G. Kruger, and Lloyd Smith, Phys. Rev. 67, 52 
(1945). 


2S. C. Sirkar and P. K. Bhattacharyya, Science and Culture 7, 89, 
(1942); Proc. Nat. Inst. Sc. India 10, 175 (1944). 


4 
| 
; 

: 


366 LETTERS TO 


On Isotopic Spin 


K. NICOLSKY 
Academy of Sciences, Moscow, U.S. S. R. 
June 2, 1945 


UPPOSE we have a dynamical system containing N 
heavy nuclear particles. A complete set of commuting 
observables for a state of system will be ra, ga, Ta Where oe 
is spin and r¢ is isotopic spin. Thus the state of a dynamical 
nuclear system is determined by the wave functions 
(rit + -tN; o1++on; Tie*+7N). For the two-valued spin 
variables oq and ry, we have q-2** functions 


Vrouw: sons tie ryt tN), (1) 


where g is the number of the independent stationary states 
determined by the Cartesian forces. 

The properties of the stationary states of the dynamical 
system with the N particles can be expressed by the 
irreducible representations of the symmetrical group of N 
variables. 

The 22" -q dimensional representation D4 with the basis 
(1) splits under the transformations of the symmetrical 
group into irreducibly invariant subspaces Dg, 


D4=D,XD.XD,=24 Da, (2) 
each defined by Partitio Numerorum 
N=K,+Ket+ sed +Ka, 
K,2=K.2=:-:-=Ka, (3) 


where K, are the whole positive numbers. 
The reduction of the D, (and D,) is determined by the 
following theorem: 


Dt== Ds", D'=z Dr, (4) 
am Tf? 


where D*”, D*" are the irreducible representations and 
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—sSm=+s, =(N/2), (N/2)—1, ---4$ or 0, 
—T=rS+T, T=(N/2), (N/2)—1, ---4o0r0, (5) 


the number of the rows and columns of the matrices 
(D*™) ag, (D™") pe; i.e., the number of independent states 
of the basis of the irreducible representation is: 


Ar, m= ((h/2)-s)-(¥/2)-s-1) }- (6) 


The Partitio Numerorum (D*™) and (D*™") are to be! 


N, (N—1)+1, (N—2)+2, 

»++((N/2)+|s|)—((N/2)—|s!). (7) 
The reduction of the D4 is thus determined by the prop- 
erties of the irreducible representation of the symmetrical 
group D® (i.e., by the representation of the Cartesian 
co-ordinates) and by the Pauli principle. 

The formulation of this problem for the p-valued spin 
and »-valued isotopic spin, i.e., for the many-valued electric 
charges, is very interesting. The physical properties of the 
dynamical quantum system are determined by the charac- 
ters X, of the irreducible representations D, of the sym- 
metrical group. 

The study of the possible types of the nuclear interaction 
is made especially easy by consideration of the possible 
irreducible representation of the symmetric group. The 
zero-approximation is determined by Hamiltonian Ho and 
the representation D,,. The interaction may be considered 
as introducing a splitting of the stationary states of the 
zero-approximation (the interchange effects for a spin and 
isotopic spin-variables). The energetic relations can be 
determined by the general Heitler- Dirac? formula by means 
of the characters X,. 

1 J. v. Neumann and E. Wigner, Zeits. f. Physik 47, 203; 49,73 (1928). 

*W. Heitler, Zeits. {. Physik 51, 805 (1928); Dirac, Proc. Roy. Soc. 123, 


714 (1929); H. Weyl, Gnippentheorie and Quantenmechanik (S. Hirzel, 
Leipzig, 1931), Vol. 2, Chap. V. 
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